{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Phonon Dispersion\n",
    "# 1. Introduction\n",
    "- In this example, we use harmonic lattice dynamics to calculate the phonon dispersion of diamond silicon."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Importing Relevant Functions\n",
    "- The inputs/outputs for GPUMD are processed using the [Atomic Simulation Environment (ASE)](https://wiki.fysik.dtu.dk/ase/) and the [thermo](https://github.com/AlexGabourie/thermo) package."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from pylab import *\n",
    "from ase.lattice.cubic import Diamond\n",
    "from ase.build import bulk\n",
    "from thermo.gpumd.preproc import add_basis, repeat\n",
    "from thermo.gpumd.io import create_basis, create_kpoints, ase_atoms_to_gpumd\n",
    "from thermo.gpumd.data import load_omega2"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 2. Preparting the Inputs\n",
    "- The structure as specified is 64-atom diamond silicon at zero temperature and zero pressure. \n",
    "- Weuse the minimal Tersoff potential [[Fan 2020]](https://doi.org/10.1088/1361-648X/ab5c5f).\n",
    "\n",
    "## Generate the  [xyz.in](https://gpumd.zheyongfan.org/index.php/The_xyz.in_input_file) file:"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Create Si Unit Cell & Add Basis"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Atoms(symbols='Si2', pbc=True, cell=[[0.0, 2.717, 2.717], [2.717, 0.0, 2.717], [2.717, 2.717, 0.0]])"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "a=5.434\n",
    "Si_UC = bulk('Si', 'diamond', a=a)\n",
    "add_basis(Si_UC)\n",
    "Si_UC"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### Transform Si to Cubic Supercell"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "Atoms(symbols='Si64', pbc=True, cell=[10.868, 10.868, 10.868])"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Create 8 atom diamond structure\n",
    "Si = repeat(Si_UC, [2,2,1])\n",
    "Si.set_cell([a, a, a])\n",
    "Si.wrap()\n",
    "\n",
    "# Complete full supercell\n",
    "Si = repeat(Si, [2,2,2])\n",
    "Si"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Write [xyz.in](https://gpumd.zheyongfan.org/index.php/The_xyz.in_input_file) File"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "ase_atoms_to_gpumd(Si, M=4, cutoff=3)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Write [basis.in](https://gpumd.zheyongfan.org/index.php/The_basis.in_input_file) File\n",
    "- The [basis.in](https://gpumd.zheyongfan.org/index.php/The_basis.in_input_file) file reads:\n",
    "```\n",
    "2\n",
    "0 28\n",
    "4 28\n",
    "0\n",
    "0\n",
    "0\n",
    "0\n",
    "1\n",
    "1\n",
    "1\n",
    "1\n",
    "...\n",
    "```\n",
    "- Here the primitive cell is chosen as the unit cell. There are only two basis atoms in the unit cell, as indicated by the number 2 in the first line.\n",
    "\n",
    "- The next two lines list the indices (0 and 4) and masses (both are 28 amu) for the two basis atoms.\n",
    "\n",
    "- The next lines map all the atoms (including the basis atoms) in the super cell to the basis atoms: atoms equivalent to atom 0 have a label 0, and atoms equivalent to atom 1 have a label 1.\n",
    "\n",
    "**Note:** The [basis.in](https://gpumd.zheyongfan.org/index.php/The_basis.in_input_file) file generated by this Jupyter notebook may look different, but the same concepts apply and the results will be the same."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "create_basis(Si)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Write [kpoints.in](https://gpumd.zheyongfan.org/index.php/The_kpoints.in_input_file) File\n",
    "- The $k$ vectors are defined in the reciprocal space with respect to the unit cell chosen in the [basis.in](https://gpumd.zheyongfan.org/index.php/The_basis.in_input_file) file.\n",
    "- We use the $\\Gamma-X-K-\\Gamma-L$ path, with 400 $k$ points in total."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "linear_path, sym_points, labels = create_kpoints(Si_UC, path='GXKGL',npoints=400)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## The <code>phonon.in</code> file:\n",
    "The <code>phonon.in</code> input file is given below:<br>\n",
    "```\n",
    "potential       potentials/tersoff/Si_Fan_2019.txt 0\n",
    "cutoff          5.0 # in units of A\n",
    "delta           0.005 # displacement in units of A\n",
    "```\n",
    "\n",
    "- The first line with the [potential](https://gpumd.zheyongfan.org/index.php/The_potential_keyword) keyword states that the potential to be used is specified in the file [Si_Fan_2019.txt](https://github.com/brucefan1983/GPUMD/blob/master/potentials/tersoff/Si_Fan_2019.txt).\n",
    "\n",
    "- The second line with the [cutoff](https://gpumd.zheyongfan.org/index.php/The_cutoff_keyword) keyword tells that the force constants will be calculated with a cutoff of 5.0 $\\mathring A$. Here, the point is that first and second nearest neighbors need to be included.\n",
    "\n",
    "- The third line with the [delta](https://gpumd.zheyongfan.org/index.php/The_delta_keyword) keyword tells that a displacement of 0.005 $\\mathring A$ will be used in the finite-displacement method."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 3. Results and Discussion"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Figure Properties"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "aw = 2\n",
    "fs = 24\n",
    "font = {'size'   : fs}\n",
    "matplotlib.rc('font', **font)\n",
    "matplotlib.rc('axes' , linewidth=aw)\n",
    "\n",
    "def set_fig_properties(ax_list):\n",
    "    tl = 8\n",
    "    tw = 2\n",
    "    tlm = 4\n",
    "    \n",
    "    for ax in ax_list:\n",
    "        ax.tick_params(which='major', length=tl, width=tw)\n",
    "        ax.tick_params(which='minor', length=tlm, width=tw)\n",
    "        ax.tick_params(which='both', axis='both', direction='in', right=True, top=True)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Plot Phonon Dispersion\n",
    "- The [omega2.out](https://gpumd.zheyongfan.org/index.php/The_omega2.out_output_file) output file is loaded and processed to create the following figure. The previously defined kpoints are used for the $x$-axis."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "nu = load_omega2()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "scrolled": false
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAoMAAAJJCAYAAADRHT1cAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nOzdd3hb5fXA8e/VsGVb3nuv7J04gwRCwg4zZZVNyixdtHRCW35QWkahtJTSCRRKGaUtZY+EQAJk7+UR24lXvKc8Za37+0O2bAc7sR3LGj6f59GjXElXemNdXR294xxFVVWEEEIIIcTEpPF0A4QQQgghhOdIMCiEEEIIMYFJMCiEEEIIMYFJMCiEEEIIMYFJMCiEEEIIMYHpPN0Ab6UoiiyzFkIIIYRPUVVVGek+0jMohBBCCDGBSc/gSUgeRjHWFMX5o02OLf8i76vwFDn2BPQdB6MhPYNCCCGEEBOYBINCCCGEEBOYBINCCCGEEBOYBINCCCGEEBOYBINCCCGEEBOYBINCCCGEEBOYpJYRYpytWLHC000QbiDvq/AUOfbEqZKeQSGEEEKICUyCQSGEEEKICUyCQSGEEEKICUyCQSGEEEKICUyCQSGEEEKICUyCQSGEEEKICUyCQSGEEEKICUyCQSGEEEKICUyCQSGEEEKICUyCQSGEEEKICUyCQSGEEEKICUyCQSGEEEKICUyCQSGEEEKICUyCQSGEEEKICUyCQSGEEEKICczrg0FFUUIVRblMUZRfKoryoaIoDYqiqD2XacN8jiBFUb6nKMomRVHqFUUxK4pS1vN8P3D3/0EIIYQQwlvpPN2AYTgHeHO0OyuKMgN4F8jquckKdAJpPZfzgCdPsY1CCCGEED7JF4JBgDpgF7ATqAT+NpydFEVJBTYAccBm4D5gs6qqDkVRQoAFwBVuabEQQgghhA/whWDwXVVV3+rdUBQlYwT7/gVnILgRuEBVVUvvHaqqdgBf9FyEEEIIISYkr58zqKqqfTT7KYoyF7ioZ/Ob/QNBIYQQQgjh5PXB4Cm4oed6n6qq+R5tiRBCCCGEl/LnYHBpz/VeRVEiFEX5jaIoJYqidCuKUqMoyv8URTnDoy0UQgghhPAwfw4GJ/f79y7gB0AyzpXE8cDlwOcnSy2jKMqQFyGEEEKI8eKumMSfg8GInus1OFPIfAsIU1U1EsgE3gcU4AlFUVZ4polCCCGEEJ7lz8Ggpt/1r1VV/ZOqqmYAVVVLgauACpwB4U+GehJVVYe8CCGEEEKMF3fFJP4cDLb3+/fvj7+zJzD8c8/mSkVRtOPSKiGEEEIIL+LPwWBVz3WTqqoNQzzmcM91EBDt/iYJIYQQQngXfw4GD43w8TLuK4QQQogJx5+DwfU911GKosQM8ZhpPddtQKP7mySEEEII4V38ORh8k755g987/k5FUQzAXT2ba1VVdYxXw4QQQgghvIVPBIOKosT0XoDIfndF9L9PURTX/0dV1Ubg0Z7NHyuK8s2eABBFUdKB/wCpgAX41fj8T4QQQgghvIvO0w0Ypvohbt963HYmUNpv+1FgOnAj8EfgKUVR2ukLKC3AGlVV9w/1whf//nP0Wg06nQaDVoNepyHUoCMhPIjUqGCyY0LIjAkh1hiITucTsbUQQggPM1vtNHV0U99moa6ti2PNZiqbu6hr68bUZcVstWO1O7A6VGw2ByqgKKCgoNMqaDUKeo1CoF5L7aTL0DgsPPRuHrGhASSEBRJtNBBjDCAmNJDo4AC0Wvl+EkPzlWBwVFRn4p2bFEV5D7gTmAcYgXLgE+A3qqrmneg5cqvbRvSaGgUCdVpCDTpiQgPJiA5mUXoUS7KimJYQikYjH0ghhPBHqqrS1GGh2mSmsqWL0oYODte0Ud7USVOHhZYuK+3dNqw9wd2YiZkKwN83l5zwYVqNQqBOQ3CAFmOgjohgPVEhASSGB5ERE8K0+FBmJYcTGRIwlq0TPsAngkFVVU+pzoqqqq8Dr49Rc07IoUKX1U6X1U5dWzd5Va18cLDGdb9OoxAWpCcjOpgF6ZGsmBLL0sxo6VUUQggfoKoq9W3dFNe3c7S+g6P1HRyubeVIXQcN7d3YHN6bmMLuUOm02Om02Glotwy5bFIBAnQajIE6oo0BpEYGMzUhlLmpESzJjCIiWIJFf+MTwaAnhQfpcThUHGrvxfmB6v33SNkczl+OTR0W9pS38NwXzl9yQXotyZFBzE+N4IKZ8aycEicBohBCeFCr2Up+VSt51a3kVbVyuKaN4vp2Oi12t7yeooBWcQ4BazUKGkVBo4BGo6AAqurMgab2fP/0fi+Zuy2AgqLRjkmPowp02xx02yw0dlgorG3nk4I61/0aBUICdcQYA0mPDmZmUhinZ8ewKCMKvXxv+SRFyqoNTlEUFThpiRdVVemw2DnW1Mm2kkYOVpg42tBBTauZ1i4rnVY7o/0TRwbrmRIfyplTYrlmUQoxRsPonkh4lZUrVwKwceNGj7ZDjC15X31bS6eFvRUtHKgwkVdtIq+6lYqmrjF9DQWICw1kcpyRGUlhzEwKZ1K8kbSoYEIN+lE/7/HHXqfFRlWLmaqWLqpauihr6KC4rp2i+naONXe5tffSoNcQawwkKzaEOSkRLJ8cQ05apMxZHAeK4hxEHc1oqgSDQxhuMHgyqqpS19ZNbmULW440sbO0iZKGDlrNthE/V2igjhlJYayalcBVOSmndPIQniNBg3+S99V32OwODte2sbe8pefSzNGGjjF/neiQAJZNiuG0rCjmpkQwOd5IoG7sK5+O5NhzOFSONXdRWNvG4do2cqtM7ClroabVfNJ9FUZfnSHUoCMtKpg5KeGsnBLHiikxGAJkcHIsSTDoBmMVDA7F1GVlb3kze8qa+aywnkOVJuwjfKnIYD1Ls6NZszSDJVlSTc9XSNDgn+R99V4Wm4MDx1rYdrSRbUeb2FPePOqhXp0G7I7BgyKtRuG0rCjOmx7PGZNjyI41ur6g3Wksjr2qli72lDezu8x5OVhpGnJUS6eBtKgQwoL0dFps1LZ202q2jngULDhAS1JEELOSwlgxJY7zZ8YTEigB4mhJMOgG7g4Gj2e22tlZ2sSmogY+O1xHQW37yXfqR69VmJ4YxlfmJXP94lT5xeXFJGjwT/K+eo/jg79dZU2YraOrKzAp1khChIGq5k6ONnR+6X6dRmHl1DgunJXAOdPjPLK4wh3HXkunhS1HGvmiqJ7PCxuobBl8yDw4QMsFMxNYPS+J5IggNhc3srusiaK6dipbumg320bUmxgSoCUjOoSFmVFcPDuBhemRkoVjmCQYdIPxDgaPV23qYn1eLevyatl6pHFEczwUICvWyJU5ydy6LEMCQy8jQYN/kvfVc1RVpaiunc8L6/mssJ5dpc10WUfe8xcSoGVeWgQL0iKZlxpBVUsXL24p5Uj9l4eQp8aHcvXCFC6fn0y0MXAs/huj5u5jT1VV8qvbWJdXw9rcWvKrWwd9XFK4gWsWpfHVRSkkhgcB4HA42H/MxMd5tewua+ZIfTtNHZZhL8BUgGhjAFMTQjljUiyr5yaRFBk0Rv8z/yLBoBt4Ohjsz9RlZX1eLW/tq2RzccOIVjErQEZMCJfPT+b25ZkES2DocRI0+Cd5X8eXqdPKpuIGPiusY0NBPfXt3SN+jozoYBakRbIgPZIFaZFMTQhFo8C6vFqeXHeYwuNGaLQahUvmJHLL6ZnMTQkflyHg4RjvY6+8sZN3D1Tx5t5Kiuu+PIqlUeDsafGsWZbOGZNiBv07Fde2sS6vlu0lTRTWtlHX1o19mF9ueq1CckQQ81IjOG9mPOdOiydQP/ZzMX2NBINu4E3BYH91bWbe21/NG3uOkVs1+K+zofT2GN6yLJ3rl6RJ17uHSNDgn+R9dS+b3cGOkibe3l/FliMNHGvqGtHwo0GvYU5KBDk9gd/8tAhijuvR23qkkcc+KmB/RcuA242BOq5bnMrXTs8kOcL7eqU8deypqkpuVStv7a3kf3sraeqwfOkx0xJCufWMTFbPSzrp4pnCmjbeO1DFliONFNe309JpHXZbwgw6JsWFsmxSNJfNTWJKfOiI/z++ToJBN/DWYLC/A8daeG1HOW/vqxrxZGidRiEnPZLvnjOZZZNi3NRCMRgJGvyTvK9jw2p3UGMy96x4bWXLkUZyq1qpauka0ahIVmwI81OdQd+81AimJoSiHyK9SW2rmYffz+ed/VUDbg8J0HLb8ixuOyOT8CDvzd7gDcdet83OutxaXttRzpYjX85mHWMM5JbTM7h5afqwM2HYbA4+K6pnXW4te8qbqWjuHPbcT51GISHcwJyUcM6dHs8FM+MJCfTe93AsSDDoBr4QDPZq77bx5t5KXthUMmh6hFCDjo5u25AnUmOglnOmx/PjVVNJjgh2c2uFN5y4xdiT93VoNruDDoudjm4bzZ0WGtstNLR3O687uqlv7eZYSxeVzV1Um0YW9IEzs8LslPCeHr9I5qVEEB588i9+q93BP7aU8ruPC+no94M6QKfhptPS+ebKbI/PBxwObzv2jta389LWMv69q+JLHRXhQXpuOyOTNcsyRhVgN7Z38/7BKjYU1JNX3Up9W/ewjxdjoI6s2BBOy4zm0rmJzEr2nqH+sSDBoBv4UjDYy+FQ2XC4jue+KGHr0S//Mpsab3QmyG4eOpFqZkwwd63I5uqcFBlGdhNvO3GLseFt76uqqj1lx7ppaO+mucNKh8VGl8VOh8VOl8VGp8WOxebArqo4HCq2nupKdoeK3UHfv3vu7199yT7gsSqqCvaebbPVTofFRme3nfZuG9220a3kHUyYQcfk+FDOmBRDTnok0xPDiA0decB2qNLED/69n8O1A+vPXzo3ifsunEaSFw4HD8Xbjr1epk4rr+4o58UtJdS2DpzTGWrQcdsZmdy+PAvjKaSTUVWVvRXNvLe/mh2lTZTUdwwI7E9Eo0B8mIGZSeGcPS2Wi+ckEh7ku6X2JBh0A18MBvs7eMzEMxuKWJtb+6X7LpqVgNGg49P8OhoGmeMBEKjTcN6MeO67aJr0Fo4xbz1xi1Pjife1o9vG4do2yho7KGvs7Ll0UNfmDABHm07Fm0SHBDAvNYJzZ8Rz6ZxEjKeYbN9qd/CnDUf4w6dFA7I0TIoz8tDqmSzL9r1pM95+TrHYHLy1t5JnNhRT3jQwPU+MMYC7z5nMtYvSCBijUnbtZisfHqxhfUEthypbqWk1D3txSpBeS0ZMMIszorhkTiILM6J8pvdQgkE38PVgsFd+dSvPfFrMB4eqByQEDdBquPG0dC6cncBznx9lY2H9kL/eM6KD+fqKbK5ZKL2FY8HbT9xidNz9vpqtdme1jIpmcquctXJLGztGXe5yPCkKhAToCNRpUBQwWx20dw9dhWlOSjirZiWwamYCWbHGMWtHcV0b3//3fg4cM7luC9Jruee8yXxtWeaYBSPjzVfOKTa7g7f3VfHMhmJKjpvSlB4dzI8umMrFsxPdEnzlV7fy7v5Kth5pori+nbZhVgFTFOd8x2nxoSzNjubiOYmkR4eMefvGggSDbuAvwWCvwzVtPLH2MOvzB/YUhgfp+cH5U7h2YQpv7K3ir58dobTxy4lVwdlbeMmcRO6/ZIZHEqv6C185cYuRGev31Wp3sLO0iS3FjWwvaWR/hQmLfWQ9fQE6Z53YGGMAkSEBGAN1BAdoCQ7ovdYSqNOi0ShoFWfqFOe/+67736bVgOa42zSKgkbDgMcG6jQYA3UE6bVUm8xsKKjjo9waCmraBm2nosCi9CgumJXABTPjSYkc29EIVVV5fWcFD7yTO+BHb056JE9ePZeMGO/8ch8uXzun2B0q/9tzjN99XEiVaWAZvMWZUTx46UxmJIW5tQ1dFjsf59ewPq+WfRUmqlu6sA6z91CrKESG6MmKNTI/NYKVU2NZnBHl8frLEgy6gb8Fg712ljbx2IcF7C5rHnD7tIRQfnHZTJZkRVPR3MljH+SzPr9u0N5CBZibGsH9l8wgJz1ynFruP3ztxO2vVFWl2+ag02Kn0+Kc1+ZwzY/rmy+n1SjotRr0Wud1gE6DXqshJFA7IFXGWLyvpk4rGwvrWJ9fx8bDdSftvdBqFLJiQpgUZyQtOpiM6BDSo4JJCDcQExpIaKBu3Ie4VFXlYKWJDw/VsPZQzZA1f3UahaXZ0ayalcB5M+KJCzW4pT1dFjs/f+sQb+w55rotQKvh++dP4Y7lWWg1vjEEeCK+ek4xW+38c2sZz2woxtTVl0ZGo8ANS9L5wflTxrXjoaS+nXf2V7GpuIGi2nZauoaf2gac1VjiwwxkxYQwOyWc0zKjWZAeOW49zhIMuoG/BoPg/D+ty6vlkQ/yKTuuF3D1vCTuv2QGMcZAHA4H/9l9jD9vHLq3MD40kDvPzOKW0zNkCHmYfPXE7e3sDpWG9m6qTWZqTF1Um8zOhROdVpo7LDR1WGjptNLcaaGj20aX1T7iVavHC9BqMBp0GAN11FSUoLGZOXf5acQYA4kNDXT2yoUGEGs0EBfm3NYcF3yYrXY+ya/jzb2VfFZYh/UERcqzY0NYnBnFnJQIZiSGMTUhFIMXJNu1O1R2lzXz0aEa1ubWDFm6LECn4czJsayalcC541C67Wh9O998Zc+AHslpCaE8de08piW4t+dpPPn6OcXUZeUPnxTx4pbSAfM4I4P1/PziGVyxINkj8/YsNjsbD9ezLq+G/RUmjjV3jaqyTaBOQ1RIAInhBtKjgpkcH8rMpDDmp0USNoYpiyQYdAN/DgZ7ma12nt9UwjOfFg84wCOC9dx/3AewrLGDX7yby2eFDYNOxA3UabhodiL3XzyDKKMMIZ+Ir5+4PanLYqesqYPShg5KGjopbeigtLGDiqZOakdQwcBT9Fpn7rPkiCASwgyYrQ42FzfQNsT8ueSIIM6aFsuy7BgWZUSNatWsu1jtDrYdbeTDQzWsy62lYYgKIMEBWs6a5qzdu3Jq3CmtHB2J9Xm1fO/1fQPmJl65IIVffWUWQQGeD6DHkr+cU4rr2njwnTw2FTcMuH355Bge/sps0qI9v5ixs9vGZ4X1fFFUz8HKVsqbOmntso4oAXp/GgUMei2hBh3hQXqiggOIDTOQFG4gJTKI1KhgkiKCiA81EBZ04p5+CQbdYCIEg72qWrp45IN83jtQPeD25ZNjeOTy2aRG9X0ALTYHv/+kiJe3lQ3o1u+l4Jz8ff8lM1iYEeXupvskfzlxu5PF5qC4rp2CmlYKatrIr26lqLadmlbzyXcegQCthqCeuXMGvRZN77y5nrlwWo2CQ1Wx2h1Y7SoWmwOr3UG3zUFHt21ENcOHa25KOOfNiOec6fFMSwj1qpWMZqudL4oa+OhQDevzawc9B4BzLvK50+NZNSuB5ZNjxrX3UlVVnt9UwsMf5LsW1wToNDx02UyuWZTqVX/PseJP5xRVVVmbW8sv38sb0MNs0Gv4/nlTuPX0THQenpt3PFVVKahu47OiOvaUtVDS0EFNq5n2btuYL/BS6Ju3q9UoKAo95ys4+ItVve2RYHCsTKRgsNdnhfX89H8HB3wAg/RafnbxdG5Ykvalk+i63BqeWHuYokFqUwLEhwXyzZWTuOk0KX3Xnz+duMdCl8XOoSoT+ytaOFRpoqCmjeK69lEFWlEhASSEGUgMN5AQbiA2NJCokAAigwOICgkgIlhPZHAARoOOYL32lL5UeucctplttHfbuH7NrTh0Bu598GEa2izUt3fT0NZNfXs39W3d1LaaaR6ivFZGdDBfmZ/MV+Yle91ihvZuGxsP1/HhoRo2FtQNmcMtxhjA+TMTuHBWAqdlRQ9Z7cOdbHYHD76by8vbyl23pUYF8ecbcpiVHD7u7Rkv/nhO6ei28duPC3lhc8mA6RwL0iJ48qvzyPSyz8lQKpo62VRUz4FjJo7UO4PE5k4Lnd127GMcX5T9+hJAgsExNRGDQXB+AJ9cV8gLW0oG/KI5a2osv75qzqCTvCuaO3nonVw2HK4f9AvcoNewel4yP794+rDLEPkzfzxxD5fdoXKkvp195S3sO9bCvvIWDte2DXt4V6tRSIkMIiM6hMyYEDKig8mICSE9OoTEcINH588N533ttNioajFT1dJFXVs3qqoyOT6UuSneVQmhpdPC+vw6PjpUw+dF9ViGSDuVHBHEBTMTWDUrgZz0SI8uxmgzW/nWq3v5vLDedVtOeiR/uynHJ6qInAp/Pqfsr2jhJ28cGDDvM0iv5acXT+fGQTopfEmXxcbBqlbyq1qdvYkmM/XtZpo7rLSarX1J4R3qsIahJRh0g4kaDPbaW97Mj/97YECvX2SwnkevmM2qWYmD7mOxOXj6kyL+OdQQsgKLMqJ48LIZzEj031/pJ+PPJ+7jWe0ODlWa2FHSxPaSJnaWNg07v1dqVBDTEsKYnhDKtMQwpsSHkhYV7LW54Hz9fa02dbEut5a1uTVsL2kaMkDPjAlh1SxnD+BsLynn1dDezc3P7yCvutV126Vzk3jiqjlescDG3Xz92DsZq93BXz87wlPrByYKP3NKLE9cNYf4MPesRPcmqqrS3m2jvs050tDQbqHLYqPb7sBic2CzO7hzxaTex0owOFYmejAIzvlBj390mL9vLhlw+3WLU3ng0pknPMmuPVTD42sLOFI/eFqJ1MggvnfuZK7MSR3TNvsCfz5xm6129le0sKOkiR2lTewua/5SbdLBTIozMjclgrmp4cxIDGNKQihhPtaL7Ivva3FdO2tza1iXW8P+fomYjzctIZQLZyWyalYCU+KNXhEA9qpq6eLG57dztN+55u5zJnPPuZO9qp3u5IvH3mgcqjTx/X/vo7C2r5MiKiSA3351LiunxnmwZd5BFpC4gQSDfTYXN/DD/+ynul9y0GkJofzxhgVkn6Q6QEl9Ow+8m8umooZB03iEBGq5ZmEqP75gKoaA8Vll6Gn+dOJ2OFTyqlvZVNzApqIGdpQ2DTmk2CvGGMi81Ajmp0UwLzWC2SnhPhf4DcYX3ldVVTlwzMTaXGcKmKF+rAHMT4twDgHPTPC6eYy9Shs6uOG57a55zhoFHrtiDl9dNLF+ZPrCsTdWzFY7v/24kGe/ODpgKtPXz8zihxdM9chcVW8hwaAbSDA4kKnLys/ePDhgxXFwgJZHLp/NV+Ynn3T/TouN36w9zOu7Kujo/nJPkVaBMybH8tDqmV5b6mes+PqJu8Zk5ouier4oamBzcQONQ9S37pUcEUROeiQzk8OYEhdKWJCOTosdm13F5lCxOxzYHWBzOIPIgJ7EzgE6jevfxkBn2oWwIL3XDvt56/vabbOzo6SJ9Xm1rMurHfCjrr/eJNDnz0zg/BnxXj/0VlDTyk3P76C+zZnSRq9V+P2187lo9uDTWPyZtx577rT1SCPf/dde6tr6UhrNT4vgD9fNH/MKNr5CgkE3kGDwy1RV5dUd5fzi3bwBvT/XL0njwUtnDnsu1393VfD7T4qoaB48MW12rJGfrJrK+TMTxqTd3sbXTtydFhvbjzbxeVE9m4oahlw93ssYqCMk0BmwdXbbabeMbXqFAJ2GMIOeiGA9scZAEsKdCZ3jQw3EhxlICA8kLtS5ong8U1B40/ta12ZmY0E9nxTUsqmoYcgVwEF6LSumxHLBrHjOnhpPeLBv9NAW1LRy3d+2uVZnG/Qa/nJjzoQdKvSmY288NbZ384P/7Gfj4b5FQxHBev5w3XyWT471YMs8Q4JBN5BgcGi5VSa+/ereAYXGF6ZH8qcbF4yopFRelYkH381jZ2nToMFCeJCeNUvTufvsyei8dNHAaHj7iVtVVQpr29l4uI6Nh+vZVdZ0wqoY3kqrUUiOCCI9Opi0qGDXdVqUcyXyWCce9uT76nA4S8B9WlDHhsN1HDjB/L+IYD3nTIvngpnxLJ8c63MJmIvr2rn2b1tpaHf2SIcG6nj+a4tYnDlx85p6+znFnRwOlWe/OMoTaw+7FpdoFPjRBdO4a0XWhJk3ChIMuoUEgyfW3m3j3jcODBg2Tggz8NebcpibGjGi5zJ1WnjkwwLe3leJ2frl+WY6jcLKqbHcd9H0k85R9AXeeOJuNVvZXNTAxsJ6NhTUDRh6GQvBAVqMgTpCDTqMBj3GQC16rQadRhmQQFUFrDYHFrszubPF5ry0ddto7bJi6rKOSWCqKM7h60lxRibFGp3XPZfRlkgb7/e1qcPC5uIGPi+sZ8Ph+iErgACkRQVz9rQ4zp8Rz+LMKK9L2jtcpQ0dfPWvW13HZ2igjpdvXzLic46/8cZzynjbXdbMN1/ZTW1r3+fgotkJPH7V3HGreuNpEgy6gQSDJ6eqzl9kj31Y4FocEqDT8Ojls7kyJ2XEz+dwOPjntnL+tLF4wAe6v6RwA2uWZXDb6Zk+21voDSduVXUu/Piop5RYUV3bqOr0ajUKqZFBpEWHkBYV5BqqjQ0LJC7UOVwbGawfs+BDVVXMVgemLmeN4bqeZM61JjO1bWZqW7upazVTZTK75pKNVIwxgOzjAsTsWCOJ4YYT9jK4+301W+3sKm3mi+J6Nhc3kFvVOuTwu06jsDAjkrOnxXH2tHiyY0N8vofkWHMn1/x1m2uxSHCAln/etpic9InbI9jLG84p3qCuzcy3XtnDztJm121T40N5bs3CAZW0/JUEg24gweDwfV5Yz3de2zsgt+B3zp7E98+bMuovoF2lTfzyvbwh0130Tna/78JpzEjyrZyFnjpxmzqtfHiomrf3VbKvwjSiguvBAVqmxIcyLSGUqQmhZMUayYh21sz01tV7XRY7Fc2dlDV2UtbYQXlTJ+VNznrG5U2dIw5+gwO0ZMcayY4NcQWI2XFGMqJDCNBpxvx9tTtU8qtb2VzcwKbiBnaUNNF9gpXaUSEBrJway9nT4lg+OZbwIN+Y/zccdW1mrvrzVsqbOgFnLfQXb1nM0uxoD7fMO0gw2Mdic/Dw+3n8Y2uZ67YYYwB/u3khC9IiPdgy95Ng0A0kGByZssYO7nxpN4dr+7LEr56XxONXzSFQN/o5SQ3tZh55v4CPcmuGzFcXHRLAxXMSufucScQYvXsFJIzfidvhUNlytJF/7yxn29GmYQ/9BgVoyUmLYF5qJLNTwpmWEEpqZDAaD1aXGGvdNjulDZ0U16G/qx8AACAASURBVLU7L/XO66P17ScMuAaj1SikRQVTU7gfvbmJH33jFlKigkmNDCI5InjYc/K6bXYOHHMm6N5R4szR2N49dIJurUZhXmoEZ0yKYcXUWOamRHi0Aoi7tHfbuOavW8mtciaUDtBqeG7NQs6cMvEWCAxFgsEv++/uY/z0fwex2HuyFOg0/ObquVw2N8nDLXMfCQbdQILBkWvvtvHtV/cMWNm1OCOKv96UQ2TI6OZh9ffBgWr+sKGIguq2IUvzpEQGcVVOCneemUWwl+YtdNeJW1VVdpY28cbuY2w92sSx5pP3fik4/2anZUWzNDuauakRZEaH+FXgNxJ2h0plcxfF9W2uQPFIfQfFde2DVtUZjhhjAMmRwSRHGIg1BhJjDCQ2NJDIkABUVeVQZSs7SpvYV9Fy0hyNWbEhLJ8UwxmTY1mSFeUX+RlPxGp3cOuLO/miqAFwBsB/uTGH82bEe7hl3kWCwcHtLG3izpd2DagJfs+5U7j7nEk+P21iMBIMuoEEg6Njszt44J1cXtneVyg+MyaEl25dPGZzNlo6LTy5rpC391XSOkRpMwVIiw7molmJ3LY8w6t6DMfqxK2qKnvKm/nvrmNsPdpIRXPXsGr8JoQFsmxSDJfOSWJRZtSEmVx9KlRVpbHDwpGe4PBIfW+g2E5lS9eYps7pLyHMwOLMKM6YFMPpk2NIjghyzwt5IVVV+cF/9vO/PZWu2x67YjbXLk7zYKu8kwSDQytr7ODWF3cOSLB+3eJUfrl6ls8upBqKBINuIMHg6KmqynNflPDIh/muL8n4sEBeunUJUxNCx/S1th9t5A+fFrG95MTpT2JCAjhzSiy3L8/0+BzD0Z64LTY7n+TX8eGhavaUt1BtMg8r+IsK1rMgPZLL5yezcmocIRL8jakui52Shg5u/s6PsRkiueDy6zjW0smx5i4qm7sG1FI9mcyYEBZnRLEoM4olmVGkRAb5ZQ/GcDyxtoA/bjji2v7uOZO557wpHmyR95Jg8MRMXVa+9coeNhU3uG47d3o8f7huvs+lVjoRCQbdQILBU/fhwWq++/o+19BXeJCev39tETnpYz+J1+Fw8N7Bap7/ooSDlaYTDo/qtQpZMUZWTo3h2kVpZI5zuprhnLhVVaWop27s1iONFNS00tQxvGHKiGA9s5PDuWhWIpfOTcTo50OJ3mKw99XuUKltNXOsuYtqUxcN7ZaeIvPdtHRaUFVIjQpmUUYUizIjR5Sn05/9a0c59/7voGv7moWpPHbl7AkbGJ+MBIMnZ7U7+MkbBwb0NOekR/L8moWjTiflbSQYdAMJBsfGliMN3PnSbtdE+CC9lr/clMMKN07+ttkcvLH3GP/eVcHBY62uCcRD0WkUEsINzEoK48wpcayYGkNyhPvSEBx/4jZ1Wdh+1Llo4MAxE0cb2mnqsAx7tWt4kJ7ZyWGsmpXIV+YlSfDnIfKFPDZ2lDRxw3PbXD39Z02N5dmbF/rdkN5YkmNveFRV5dcfHeYvn/X1OE+KM/LybUtICPf9H2ISDLqBBINj5+AxE2te2EFTTw1bvVbh6Wvnc+E41RDdUtzAC5tL2F7SNOQcw+NpFAg16EgICyIxwkBKRBCZMSFMSQglKSKIuNBAQgK0aDRDf0HZbA7q2s3UmLqpa+umpL6dksZO3ln/BfaAEIwxSbSabcMa6u2lKBBnDGR2SjgXzEzgotmJMuzrJeQL+dRVNHWy+o+bXeeK6Ylh/PeupXKMn4QceyPz900l/PL9PNc0prSoYF69Y4nP1zSWYNANJBgcW0fq27n5+R2uhLFajcLvr53HJXPGd5l/S6eF13dWsC6vhvzqtiHT1QzXYJ+4sTxiggO0pEQGMT8tkkvnJLI0O8Yv04f4A/lCPjUd3Tau/PMWCmqc6amiQwJ4+9un+/wX9HiQY2/k3t1fxff/vc/VA50UbuDVO04jIybEwy0bPQkG3UCCwbFXberihue2c7RnVZdGgd9dM4/V85I91iZTp4WPcmvYeLieQ5Umatu6T5rewx20GoUwg46UyCBmJ0ewcmosyyfHEOSl6XHEl8kX8ug5HCp3vbybdXm1gHP04LU7TmNhhlQXGQ459kbnk/xavvHyHtdUorjQQF65fQmT48d2oeN4kWDQDSQYdI+6VjPXP7ed4rp2wBkQ/ubquVyxYOTl69zFbLGxs6yZ7UebyK9upcZkpqnTQqvZitnqwOFQR9T7p1GcwZ5BryXMoKe+4ghaSztrLl9FToZz1WiYH1WLmKjkC3n0nlpfyFPri1zbj185h68uSvVgi3yLHHuj90VRPXe8tAuz1RkQRoUE8PJtS5iRFObhlo2cBINuIMGg+9S3dXPjc9td1UoUBZ64ai5XjaKesSe1m23Ut5nptjnQahQ0GgWtoqDVOC/RIQEE6r+ctkBO3P5J3tfR+aywnq+9sMM1f+vW0zP5v0tneLZRPkaOvVOz/Wgjt764k46eaUPhQXpeunUxc1MjPNyykTmVYFCWZ4lxFxsayKt3LGFaT85BVYUf/3c/7x2o8nDLRsZo0JEZa2RaYhiT40PJjjWSERNCapSzZu9ggaAQok9lSxff+9deVyC4NCuan140zbONEhPOkqxo/nn7EkINzmk5pi4rNz6/nUOVJg+3bPxIMCg8ItoYyGt3nMb0RGdXvEOF7/1rH+t75gwJIfybxebgW6/scZUKiwsN5Onr5ksKGeERC9Iiee2O04gMdk7ZaTPbuOn57RzuWdDk7+RTJzwmMiSAf962mOxY5+otm0Plm6/uYVNRw0n2FEL4ukc+yGdfRQvgnFP7zPULiA0N9HCrxEQ2KzmcV24/jfCeOdzNnVZueG6ba467P5NgUHhUjDGQV24/jbSeusUWm4M7XtrFztImD7dMCOEu7+6v4sUtpa7te1dNY3GmrBwWnjcjKYyXbl1MaE9uy4Z2Czc8t42yxo6T7OnbJBgUHpcQbuCV25eQ1JMBvstq59YXd1JQ0+rhlgkhxlp5Yyf39Ss1t2pmArcvz/Rgi4QYaG5qBC/euojgnrrFta3dXP/sdo41d3q4Ze7j9cGgoiihiqJcpijKLxVF+VBRlAZFUdSey4hnGiuK8lS//Te6ocliFFKjgnn59iXEGJ3DRG1mG2v+vsOvP3xCTDRWu4O7/7XXVZ4yPTqYx6+eIzWHhdfJSY/i+TWLCNQ5w6TKli6uf3Y7ta1mD7fMPbw+GATOAd4Gfg6sAqJH+0SKouQA3x6jdokxlhVr5B+3LsLY0z1f29rNzX/fQXNPaSohhG97an2ha56gTuMsSxkmtbSFl1qaHc2zNy8koGdRU3lTJ2v+vgNTl9XDLRt7vhAMAtQBHwC/AO4czRMoiqIB/oqzWtjusWuaGEszk8L52805rg/f0foObv3HTjotw6spLITwTluPNPKnjUdc2z84f6rP5XETE8+ZU2L5840L0PWUAS2oaetJUn1qpUy9jS8Eg++qqhqvqurFqqo+CHw8yuf5DpAD/AE4NFaNE2NvWXYMv7tmHr0jR3vLW/j2q3ux2ce/TJwQ4tQ1d1i45/V9rnyCy7Kj+fqZWZ5tlBDDdM70eH595RzX9o6SJu5+bS92h/8UpfD6YFBV1VMOvxVFSQF+CVQBD5xyo4TbXTwnkQcu6atC8GlBHb96P9+DLRJCjIaqqtz7vwPU9My1igzW87tr5qHRyDxB4TuuzEkZkBB9XV4tP3/roN9UKfP6YHCMPA2EAt9XVXViZJD0A187PZNvrMx2bb+4pZQXN5d4sEVCiJF6Y08la3P7ksk/cdVc4sMMHmyREKNz55nZ3NmvR/u1HRX87uNCD7Zo7Ph9MKgoyqXA5cB6VVVf93R7xMj86PypXDw70bX90Ht5bCio82CLhBDDVdnSxS/eyXVt37AkjXNnxHuwRUKcmntXTeOK+cmu7ac/LealraUea89Y8etgUFGUEOAZwAJ8a5TPMeRFuJ9Go/DkV+cyr2eiuUOFb7+6h/xqyUEohDdzOFR+/N/9tPVLI/Ozi6d7uFVCnBqNRuHXV81h5dRY120PvpM7bqVU3RWT+HUwCDwEpAFPqKrqH325E5BBr+XZmxeSHBEEQIfFzm0v7qSuzT/zPQnhD/65rYzNxY0AaBT47VfnEhyg83CrhDh1eq2GP92wYEAnxXde28vBYyYPt2z0/DYYVBRlHvBdoBR4eLTPo6rqkBcxfmJDA3nhlkWuEkFVJjPffHkPFpusMBbC2xytb+fRD/sWfN15ZjY56VJuTviP4AAdz61ZSEqks5Oiy2rntn/spKqly62v666YxC+DwZ6cgn8DtMDdqqq6990R42JKfCjP3LCA3kWIu8qa+cW7uSfeSQgxruwOlR/8Zz9mq/OH2rSEUO45b7KHWyXE2IsxBvLiLYsIMzg7Keraurn1xZ20mX0vKbVfBoPAGmARsA7YoCiKsf8F6B2r0Pa7Xeux1ophWzEllnsv7Fve/8r2cl7dXu7BFgkh+nthcwl7y51VRvRa55zfQJ2cXoV/mhQXyl9uyhmQlNoX8+L6azCY3nN9PtA2yOWGnvvP6Hfb8nFuoxilO5ZncdncJNf2A+8cYldpkwdbJIQAKGvs4DfrDru2v3P2ZGYmhXuwRUK437LsGB7rl5T6s8J6/u+dXJ+aTuavwaDwY4qi8Osr5zAjMQwAq13lrpf3UGOSBSVCeIqqqtz7xsEBw8P984QK4c+uyknh7rMnubZf3V7O85t8Jy+uXwaDqqo+qKqqMtQF+EfPQz/rd/tGDzZZjFBQgJa/3ZxDVEgAAA3t3Xzr1T1YfaxrXgh/8frOCrYe7Vs9/MRVc9Fr/fIrRohB3XPeFFbP6xu1euSDfD4rrPdgi4bPJz6piqLE9F6AyH53RfS/r2fhiJggUiKD+eP1C9D2zNXYXdbMrz8s8HCrhJh4akxmHu5XLvKOM7OYnSLDw2JiURSFx6+aw4K0gXlxj9S3e7hlJ+crwVN9v8uefrdvPe6+tPFvmvCkpdnR/OiCqa7t5zaV8NGhag+2SIiJRVVVfv7WIVdy6YzoYO45d4qHWyWEZwTqtPzlphwSw50lF9vMNm7/xy5Mnd69wthXgkEhhnTn8izOnd5X4upH/zlAaUOHB1skxMTx/sFq1uf3VV947Mo5GPSyelhMXHGhBp69eSEGvTPEKmno4Nuv7fHqFcY+EQyeaP7fcZfSYT7f13oev9K9LRfjQaNRePLquaRGOZN/tnXb+MYrezBb7R5umRD+rbnDwgNvD6w9fFpWtAdbJIR3mJUczpNXz3Ntf1HUwMMf5J9gD8/yiWBQiJMJD9bzp+tzCOiZsJ5f3SoJqYVws1+9n09jhwWAxHDDgBygQkx0F89J5O5z+hKuv7C5lNd3emdeXAkGhd+YnRLOA5fNcG2/tqOC9w/I/EEh3GHb0Ube2HPMtf3w5bMINeg92CIhvM/3zpnMhbMSXNs/f+sQO70wL64Eg8KvXL84jUv7JaS+938HONbc6cEWCeF/LDYH9791yLV90ewEzp4Wf4I9hJiYNBpnFZ7p/fLifv2fu6lo8q7vJQkGhV9RFIWHL5/lKh7eZrbx3X/t8+qJu0L4muc3lVBU50yXERKg5f8umenhFgnhvYIDdDx7cw7RPXlxmzos3PHSLjp6VuB7AwkGhd8JM+h5+rr5A/IPPv1JkYdbJYR/ONbcOeDzdM95U0joSaMhhBhcSmQwf70pB722r4bx9/+9D4fDO0rWSTAo/NKCtEi+f15frrNnNhSzrac6ghBi9H7xbh5dPSv1pyWE8rVlGZ5tkBA+YmFGFA9/ZbZre21uLU95SUeFBIPCb921Iptl2c40Fw4V7nl9H6Yu7078KYQ3W59Xy8d5fTkFH758FjopOSfEsH11USq3nJ7h2n76kyKvWOgon2Lht7Qahd9dM4/IYOcKx2qTmQfePnSSvYQQg+my2Hngnb50TdcuSiUnPcqDLRLCN/3souksnxzj2v7Bf/ZxqNLkwRZJMCj8XHyYgUevmOPafmtflVf8ChPC1/zh0yIqW7oAiAzW85NVklNQiNHQaTU8c90CMmNCADBbHdz50i7q27o91iYJBoXfWzUrgSsXpLi2f/bWQepazR5skRC+pbiujWe/OOravu/C6UT2rIwUQoxceLCeZ29eSGigDoAqk5m7Xt5Nt80zlbMkGBQTwgOXzSA5wplupqXTyo/fOICqescqLiG8maqq/PytQ1jtzs9LTnokV+WknGQvIcTJTIoz8vT18+lJfMHusmZ+/uYhj3w3STAoJoQwg57fXD0XpedDt/FwPa9s986yQEJ4k7f3VbHtqLNiglaj8KuvzELT++0lhDglZ02N474Lp7u2/7P7GH/fXDru7ZBgUEwYS7Ojue30TNf2w+/nU97oXVnghfAmpi4rv3o/z7V9y7IMVyUFIcTYuH15JlcsSHZtP/x+HhsP141rGyQYFBPKDy+YypR4IwBdVjs/keFiIYb0m7WHaWi3AJAQZuB7/XJ3CiHGhqIoPHL5bOanRQDOVGjfemXPuK4wlmBQTCgGvZYnrprrmqOx9Wgjr+2o8GyjhPBCB4618PL2Mtf2A5fOwNgz2V0IMbYMei1/vSmHpJ5qPh0WO197YQdljR3j8voSDIoJZ25qBHecmeXafuSDfKpNXR5skRDexe5Q+dmbh+jtNF8xJZZVsxI82ygh/FxcqIF/3LqY8CBnbtyGdgtr/r6Dhnb3p5yRYFBMSPecO8WV46m928ZP/3dQhouF6PHq9jIO9gxRBeg0PLR6Jooii0aEcLfJ8aE8t2YhgTpneFba2MltL+6ko9vm1teVYFBMSAa9lsevmuNaXbzhcD1v7av0bKOE8AJ1bWYeX3vYtf2tlZNIjw7xYIuEmFgWZUTx9HV9KWf2HzNx5z930WVxXw5CCQbFhLUoI4qbT0t3bf/i3TyPZoAXwhs8+kEBbWZnL0RmTAhfX5F1kj2EEGPtgpkJPLR6lmt7c3Ejt7y4w209hBIMigntx6umDUhG/cA7UrtYTFxbjjTw5t6+HvKHVs/EoNd6sEVCTFw3npbOjy6Y6tredrSJNX/fQavZOuavJcGgmNBCAnU8duVs1/YHB2v46JDULhYTj8Xm4P63+n4MXTInkeWTYz3YIiHEt86axH0X9tUB31XWzJV/2kJF09jmyJVgUEx4yyfH8tWFfeW1/u/tXNrc8MtLCG/27BdHOVLvTGNhDNRx/yUzPNwiIQTA11dk83/9Po9Fde1c/PQXfHhw7DouJBgUAvjZxTOIDQ0EoK6tmyfXFXq4RUKMn4qmTv7waZFr+/vnTSE+zODBFgkh+rv1jEx++9W5BGidYVur2cY3XtnD7f/YRUFN6yk/vwSDQgDhQfoBPSEvbS0d1+zvQniKqqr84t1czFYHADMSw7h5afpJ9hJCjLcrFqTw2p2nuea5A6zPr2XVU1+w+plNp/TcEgwK0ePSOYksnxwDOMsB/ezNg9gdkntQ+Le1ubWsz++rg/qry2eh08pXgxDeKCc9kg/uXj5gahM408+cCvnEC9FDURQeWj2LgJ5kn/uPmXi1XzkuIfxNe7eNX7yb69q+fkkaC9IiPdgiIcTJhAfrefyqubzz7dO5eHYiWs2pJ4SXYFCIfjJjQvjmymzX9uNrD1PXZvZgi4Rwn999XEi1yXl8xxgD+MkF006yhxDCW8xJieCPNyxgz8/P4883LDil55JgUIjj3LUi21Wqrs1s4+H38z3cIiHG3qFKEy9sLnFt//ziGYQH6z3YIiHEaIQH67lwduIpPYcEg0Icx6DX8st+md/f3lfFpqIGD7ZIiLFld6j87K1D9E6JPX1SNKvnJXm2UUIIj5FgUIhBnDE5hsvm9n053v/2IcxW99WFFGI8vbqjnP0VLQAEaDX8cvUsFOXU5x0JIXyTBINCDOHnl0wnNFAHQElDB89vKjnJHkJ4v7o2M49/VODa/sbKbLJijR5skRDC0yQYFGIIcaEGftivLuQfNxRTY5LFJMK3/eq9fNrMzmL3mTEhfKPfgikhxMQkwaAQJ3DDkjSmJYQC0Gmx89iHsphE+K7PCut5Z3+Va/uXq2dh0Gs92CIhhDeQYFCIE9BpNfzfpX2VSd7aV8XusiYPtkiI0WkzW7nvjQOu7dXzkjijJ8m6EGJik2BQiJNYlh3DhbMSXNsPvpOHQyqTCB/z648KqOqZ5hAZrB9Q+F4IMbFJMCjEMPz0oumuyiQHK038d/cxD7dIiOHbeqSRl7eVu7YfvGwm0cZAD7ZICOFNJBgUYhhSo4L5+plZru3H1xbQarZ6sEVCDE+nxcZP+g0Pnzs9fkDaJCGEkGBQiGH6xspsEsMNADS0W/jDJ0UebpEQJ/fkukLKmzoBCDXoePhyySkohBhIgkEhhik4QMe9F/bVbn1hcylH6ts92CIhTmx3WTN/71dy7v5LZhAfZvBgi4QQ3kjn6QYI4Usum5vEy9vK2FnajM2h8sv38njxlsWjfr6WTgsFNW2UNXZQ0dRFVUsXLV1W2sxWWrtstJqt2B0qWo2CRlHQaMCg0xIZHEBEsJ6okACiQgJIiQwmLSqY1KggkiKC0Gvld95E195t457X96H2rHVaPjmGq3NSPNsoIYRXkmBQiBFQFIUHLp3Jpc9sQlVh4+F6Pi+s58wpsSfd12Z3kFfdiilhIeawFJY++gnVbkhirdUoZEQHMy0xjOkJoUxNCGN2cjgJ4dIjNJH84p3cvuHhQB2PXjFbhoeFEIOSYFCIEZqVHM7VOSn8e5dzRfEjH+Rz+qQYtJovf9G2d9vYeLiO9Xm1bDhcj6nLChlnAdDlpmomdofKkfoOjtR38P6BatftieEGFqRFsiA9kgVpEcxMCnetkBb+5cOD1fyn34r3X10+i5TIYA+2SAjhzSQYFGIUfnD+VN7dX02X1U5BTRv/3V3BNYvSAGcw9kVRPW/urWRtbg1mq2PI5wnQaZgSbyQ71khqZDApkUFEGwMJM+gINegJNegI0GmwO1TsDhVVhU6rjeYOK82dFpo7LdS1dlPR3ElFUycVTV3UtA4eZFabzLx/sJr3D1a7XjsnLZIzJsdw+qQYZieHDxrQCt9SberivjcPurZXz0ti9bxkD7ZICOHtJBgUYhTiwwx8fUUWT613rih+Yu1hlmVH8/7BGv65tYzKlq4h9guk/eheDK0VvPr0w2THhqAb4/l9nRYbRbXtFNS0kl/dRl51KwePmeiy2gc8zmJzsPVoI1uPNvLE2sOEGXQszY7mjEnO4DAzJkSGFX2MxebgW6/soaXTmfYoOSKIh1bP8nCrhBDeToJBIUbpzjOzeG1HObWt3TS0W1j++MZBHzc1PpQLZydw7vR4ZiaFcdZZv3Le3lPzeKwFB+iYmxrB3NQI1202u4OCmjb2ljezp7yFPeXNlDV2Dtiv1WxjbW4ta3NrAUiNCuKcafGcPS2OJVlRBOqkhq23e/TDfPaUtwDOuaO/u2Ye4UF6D7dKCOHtvD4YVBQlFDgLWAQs7LmO7rl7uqqqBUPsZwAuBi4EFgNZgB6oBbYCf1ZVdaNbGy/8WnCAjhVTYl1zB/vTahTOmhrLPedOYUZSmMd72HRaDbOSw5mVHM5NS5231baa2VzcwKbiBjYXN1Db2j1gn4qmLl7cUsqLW0oJDtByxqQYzpkex1lT44iT9CRe570DVbywudS1/ZNVU1mcGeW5BgkhfIbXB4PAOcCbo9jvXeDcftvdgBVI7bl8VVGU36uq+r1Tb6KYiMxWO58X1g96n92hsj6/jmqTmfsvmcFpWdGDPs6T4sMMXLEghSsWpKCqzkUnvcHh1iONtHfbXI/ttNhZl1fLujxnr+Hs5HDOnhbHudPjmZXs+WB3osutMvGT//ZVGTl/Rjx3LM86wR5CCNHHF4JBgDpgF7ATqAT+Nox99EAR8Czwbm8PoqIo2cCjwNXAdxVFKVRV9U9uabXway9uKaWmpzctxhjAG99YxsvbyvjntjLXopHcqlau/ds2Lp6dyH0XTfPaFZ2KojApzsikOCNrlmVgsTnYWdrEJ/l1fFpQS+lxQ8oHK00crDTx+0+KiA8L5Oxp8Zw3I45l2TEY9DKcPJ6qTV3c+uJOOizOOaHp0cE8cfVcCdCFEMOmqL0ZSb2UoihaVVXt/bYzgN6U+icaJl4GbO+/b7/7FGA9cDZQoqrql35CK4qiAnj730d4RkunhTMf30Cr2dl79svVM7lpaQYA9W3dPPvFUV7aWjpgJXGgTsNdK7J585FvolFtbNy4cfwbPkpH69v5tKCOT/Lr2FnahM0x+OfCoNdwxqRYzp0ex9nT44gLnTjDyStXrgQY1/e1zWzl6r9spaCmDXCWm3vjG8uYEu+e+ajCO3ni2BPep/cHoKqqI/4l6PU9g4MFc8Pcb8sJ7lMVRXkJZzCYqShKlKqqTaNto5h4/rih2BUIZsaEcO3iNNd9saGB/PSi6dxyegaPfVjA2/uqAOi2Ofj9J0Xo5qwh5shHHmn3aGXFGsmKNXL78ixMXVY+L6zn04I6Nhyuc61cBTBbHazPr2V9vnM4eW5KOOdMj+fc6fFMTwyV3qox1GWxc9fLu12BoE6j8JcbcyQQFEKMmNcHg27U2O/fMq4lhq2qpYt/bC1zbf/ogqmDln9LDA/i99fO56bT0nnw3VwOVbYCYAuKombmdTz4Ti4/umAqIYG+9TEMD9Jz6dwkLp2bhM3uYE95C5/k1/Jxfi1H6zsGPHb/MRP7j5n47ceFJIUbOGd6POdMj2NpdrSsTj4FnRYbt/9jF1uO9J3GHr1iNqdPivFgq4QQvsrrh4mPN9xh4mE8zxPAD3GuLk5Uj/tDyDCxGMp9/zvAazsqAJibGsFb31x20h4vu0PltR3lPPZhwYCFGSmRQTxx1VyWZnvfApPRKGno4JOensGdpc3YhxhODg7QsnxyDOdMd6auiTEGjnNLx954DdV1Wmzc+uJOth3tG8z4+6xz1QAAIABJREFU0QVT+dZZk9z6usJ7yTCxgFMbJp6QwaCiKMlAAWAEfq2q6r2DPOakfxhf+9uJU1fS0MG5v/3MFeS8cvuSEfXGVLV0cc5PnqMrsm+aqqLAN1dm871zpwzaw+irTJ1WNhbWsT6/jo2H62gz2wZ9nKLAzKQwTp8UwxmTYliUEeWTi1DG4wu5oqmTr/9zN3nVra7bJBAUEgxOHMOZauOXcwbHmqIoOuAVnIFgOc6VxUIMy1PrC12B4NKs6BEPyyVFBBF3+A06YmZgnf0VTF1WVBX+uOEIm4sbefra+aRFe+eK45EKD9azel4yq+clY7X3rU7+JH/g6mRVhUOVrRyqbOWvnx0lQKdhYXqkKzicmRQ25lVafNHm4ga+/eoemvvN0bz3wmnctSLbg60SQviDCdczqCjKn4G7AAtwtqqqm4d4nAwTiwHyq1u56Okv6D0k3vjGMnLSI0f8PL2/4l9/5yN+8O/9bCpucN1nDNTx8OWz/LqWbG9Ow0/ya/kkv47d5UMPJ4NzSHleagQL0yPJyYhifloEYQbvq6rhrt6ZLoudpz4p5NnPj9L7Z9JrFR5aPYvr+i1cEhOX9AwK8PPVxGNJUZRHcAaCduCGoQJBIQbz5LpCVyB47vS4UQWC/cWHGXjp1sU8+8VRnlh7GJtDpb3bxnf/tY+dpU3cf8kMv1xk0T+n4ddXZNNqtrL9aJMr4XVxXfuAx3da7Gw50uhaLKEoMDnOyKykcGYkhTGz59rfyq512+y8vbeK360vpNpkdt0eFxrIn2/MOeXjTwghek2YYFBRlJ8B9wEqcIeqqv/1cJOED9lT3uxKlwLw/fOmjsnzajQKX1+RzdLsaO5+ba9r+PTlbeUcqmzlTzcsICkiaExey1uFGfScNyOe82bEA1BjMrPliDMw3Hakkap+gRA4h5ULa9sprG3nf3srXbenRgUxJS6UzJgQMmJCyOq5TggzoNH4Rkobs9XOvooW1ufV8vb+KurbBpYIXJYdzVPXzJNygEKIMTUhgkFFUe4BftWz+V1VVV/wZHuE73ly3WHXvy+dm8SMpLAxff45KRG8d/dyfvLfA7x/sBqAfRUtXPKHTTx97XzOmDxxUoYkhPeVyQPnoptdZc3sLm1iV1kz+dWtDDaqXNHURUVT15duD9RpiA8zEBcaSFxYIHGhBuLCAokOCSAkUEdIoA5jv4tOq6BVFBRFQatR0CigoGCxO7DYHVhtDqx2B9091xabA7PNQWfkJBwaHf/ZVYHZ5qDbaqe759psc2C22um2OjDbBl47VBWr3UFjh4VjzV2DDpnHGAP44flTuWZRquRqFEKMOb+fM6goyjeA3nJz96qq+uthvo7MGRSAc+L+Dc9tB0CrUfj4njPJijWO+vlONL9HVVWe31TCox8WuIICjQI/vWg6t52RKYEA0NFtI6+6ldxKE7lVreRWtVJU14bV7n+f1fiwQNYsy+DmpRkYfSwfpRg/MmdQgMwZHJKiKGuAP/ZsPjTcQFCIXqqq8sTavl7BqxaknFIgeDKKonD78v9n777DoyyzBg7/nknvIQkJpBF6771XxQYqiG3tva66rnVd111d3c+y6roudrGLKAKCHaT3EnoPIZVU0ttk5vn+mMnMBCmBzGSSybmvK1fyTntP0Myc9ynndKJffDj3fr6VvNJqzBqeW7KXAzmlPHtZH49cR3g2gvy8GZoUwdCkCNttNbVmDuWWkZJfRmp+OSn55RzJLyc1v7ze7tvmTinoGBnE8E6RTO0dw+guUR5VbkgI0Ty1iGRQKeU4R+a4ajr8hPsKtdZm63NmAu8DCnhJa/0310cqPM3SvbkkpxcB4Otl4I9TujbJeYd1jGDJ/WO469MtbE2znP+rzRkcyS9n9nWDPaJIszP5ehvoFRt60un70iojuaXV5JRUkVdaTW5JNbmlVRRVGCmrrqWsupZy23cTtWYzZg1ms8asNSazRmOZbvbxsnz5Wn/29VL4eBnw9/Fiy8b1KF3LBedNxt/HCz9vw0m/+/sY8PO2f1cKvA0GwgN9SGgTSIBv6072hRBNr0VMEzekALRVR611qvU5KUBH6+05p3yGxYwTexnLNLHQWnPJG6vZnWUp8HvTqCSemd670a97NlM61bUmnpi/k/lb7Rsl4sIDeO/GIfRs79x1i6JxZKpOuIv8vyegcdPEnjz/4Pi7xZzhy7fJoxPN3tK9ubZE0M/bwD0Tm764r5+3F6/M6s+TF/WgbrlgZlElM2ev5efdx5o8HiGEEJ6nRSSDWmvVwK9Uh+ckncXzlrvvtxPNkdaa15cetB1fN6ID0SHuKeehlOKOcZ1574Yhtk0EFTUm7vx0C++vPnKGZwshhBCn1yKSQSGa2rJ9uezMLAYso4J3ju90hme43uSeMcy/ZxSJEZZ2dVrDs4v38Myi3aft4CGEEEKcjiSDQpzgxFHBa4cnum1U8ETdYkL49p5RDEoMt902Z20qd326hcoakxsjE0II0VJJMijECZbvz2NHhn1U8O7xTb9W8HQig/34/PYRXNinne22X/bkcPU7637XsUIIIYQ4E0kGhXCgteY1h1HBa4YlNsvWX/4+Xrx57SDuGGefvt6eUczl/1vDodxSN0YmhBCipZFkUAgHKw7ksb2urqC3gbsnNK9RQUcGg+LJi3ry7KW9qWu9m3G8khn/W8v6lAL3BieEEKJJVBlN/LircdUlWkTRaSGawolrBa8ZmkBMMxwVPNH1I5OIDQ/gvs+3UWk0UVJVyw3vb+SlWf24dECcu8MTQgjhAuXVtcxefpg5a1Mpq65t1GvJyKAQVqsO5rMtzd5t5O4JXdwcUcNN7hnDV3eOpG2IpTNJjcnMA18m899lB6VwuhBCeJitacc5/9WV/Pe3Q41OBEGSQSGA348KXj0sgXZhzX9U0FHf+DC+vWcUXaPtvZNf/vkAj3+zE6PJ7MbIhBBCOMuCbZlc9fY6Mosqbbd1jApq1GtKMigEsO5wAVuOHgfqRgWb71rB04lvE8jXd49iZKdI221zN6dzy5xNlFYZ3RiZEEKIxvpiYxoPfZWM0WSZ8Qn19+blWf1Z+qfxjXpdSQaFAP63/LDt51lD4mkfFuDGaBonLMCHj24ZxoxB9vWCqw7mM+utdWQXV57mmUIIIZqrJTuyefLbndSt/OkeE8KSP47lisHxGAxn3Y64HkkGRau3Pb2I1YfyAfAyKO5qZnUFz4Wvt4FXZvXngcldbbftO1bKZW+uYXdWsRsjE0IIcbaS04t4aG6yLRHsFx/G3DtHkGDtSNVYkgyKVu9/yw/Zfp7Wr73T/rjcTSnFQ+d14+VZ/fG2XjXmlFRz5VvrWL4/183RCSGEaIjj5TXc+9lWaqxrvzu1DWLOzcMID/R12jkkGRSt2sGcUn7anWM7bkk7iBvqisHxfHTLMEL8LZWkymtM3PrRZr7YmObmyIQQQpyO2ax56Ktk22aREH9v5tw0jIgg5yWCIMmgaOVmr7CvFZzSM4bu7ULcGI3rjO4SxTd3jyIu3LIW0mTWPDF/Jy/+uA+zWUrPCCFEczR7xWGW78+zHb8yqz+Jkc6fvZJkULRa6YUVLEzOsh3fM7HlrxU8nW4xIXx7zyj6xIXabvvf8sM8MDeZKqPJjZEJIYQ40a7MYv79ywHb8R3jOnF+73aneca5k2RQtFrvrkrBZB0VG9kpkkGJbdwcketFh/oz946RTOoRbbvtu+1ZXPXOenJKqtwYmRBCiDpGk5lHvt5h+4wa3KENj0zt7rLzSTIoWqW80mrmbkq3HXv6qKCjID9v3rl+MNeP6GC7bXt6EdPeWE2ytS+zEEII93lr+WH2ZpcA4O9jqQ7h4+W6lE2SQdEqfbDmCNW1lp1ZfePCGNMlys0RNS1vLwP/uLQ3f5/eGy/rTuPc0mqufHsd87dmuDk6IYRovQ7mlPLGMnuVi4fP605SIzuMnIkkg6LVKaky8um6o7bjeyZ0RqnGFexsiZRS3DgqiY9vGUZ4oA8ANbVm/vTVdp7/fq9tekIIIUTTMJs1j8/faSsj0z8hnFvGdHT5eSUZFK3OJ+uOUmpt7N25bRBTXbQgt6UY3SWKhfeOrtfT+J2VKdz04UYKy2vcGJkQQrQuC5Izba1RfbwUL13RzzZ740qSDIpWpcpo4sM1qbbju8Z3bnQbH0/QITKI+feMYkpP+8aSVQfzufg/q9hytNCNkXkWrTVrD+dz72dbGfrPXxn87C/c9ckWKmpq3R2aEMLNyqpreeGHfbbj28Z2oltM05Q7k2RQtCoLtmWSX1YNQLtQfy4dEHeGZ7QeIf4+vHP9EO6baC+8nV1cxVVvr+e9VSloLdPG56qkysicNUc479WVXPvuBpbszCavtJqC8hp+3H2Mx77ZKf++QrRybyw7SF6p5fMpJtSv3nuxq3k32ZmEcDOzWfPuqhTb8S1jkvD1lushRwaD4s9TuzOoQzgPzd1OcaWRWrPmuSV72XikkH/N7Of0yveebE9WCZ+sP8rC5Ewqak5dy/G77Vn0jw/jtrGdmjA6IURzcTivjA9WH7EdP3FhT4L8mi5Fk09C0Wr8tj+Xw3nlAAT7eXP1sEQ3R9R8TeoRw5I/jqF/Qrjttp/35DD1tZX8Jn2NT6u61sSCbZnMnL2Wi/6zii82ptVLBIP9vLlhZAe+u28M0/rH2m5//vu9rD2U746QhRBu9q8f9mE02WsKXjog9gzPcC4ZGRStxjsr7aOC1wxLINTfx43RNH/xbQKZd+dInv9+L3PWpgKW+ow3f7iJ60Yk8uRFPQn0lbeQOka/MMqi+zHyhWUn3XjTPSaE60d24LKBcQRbr/hfmdWfzOMVbE0rwqzhvi+2sei+0cS3cX67KSFE87Q5tZBf9uTYjp+Z1rvJK1zIyKBoFbanF7HhiGUjhLdBcfNo12/V9wS+3gaemd6bD24aQlSwn+32T9encfF/VrP2cOseyaquNbF4RxbXv7+BzAG3Uxw3ol4i6OOlmN4/lnl3jeTHB8dy3YgOtkQQLP++s68bTNsQy79tYXkNd326RdoDCtFKaK35vx/tm0YuHRBL3/iwJo9DLutFq+C4VvCSfu2JDQ9wYzQtz6QeMfz0YDhPfruTn3ZbrmCP5Jdz7bsbmDEojicv6lkvWfR0B3JKmbspnflbMzheYbTc6HAlHxvmzx9GdODKIQm2RO9UYkL9mf2HQVz9znpqzZpdmSU8+e1OXpnVv1XWvxSiNVm2L5dNqfZSMg+f57qWc6cjyaDweOmFFXy/M9t2fPs4WaR/LiKD/XjrusF8szWTZxbtpsxaq3H+1kyW7s3l8Qt7cNWQBI8t1VNeXcuSHdl8uSmNrWknadunNf7Fqfzn/iuY1CMa77NoHTUkKYK/TevFXxfuBiz/pv3jw7lxVJKTohdCNDcms+bFH/fbjq8dlkhipHuWiEgyKDzeB2uOUNdMY0yXKHrHNv0QvKdQSnHF4HjGdo3iH4v3sGSHJckurjTyxPydfLr+KE9c2JMxXT2jvZ/JrFl3uIBvt2Xy465syk+yIzguPIBZQ+L5+sWH8a4p5fze953Tua4b0YEdGcXM22JpB/js4j30bB/KsI4RjfodhBDN04JtmezPKQUg0NeL+yZ1dVsskgwKj1ZcYWTupnTbsYwKOkdMqD9vXjuIWYNzeXrhbtIKKwDYnVXCde9vYGzXKB6c0o3BHdq4OdKzp7Vmd1YJC7Zlsmh7FrnWul+OfLwU5/dqx5VDExjTJQovg2LBc6WNOq9Simcv68P+nFJ2ZBRTa9bc89kWFt8/lnZh/o16bSFE81JlNPHvXw7Yjm8b2+mMS0pcSZJB4dE+23jUVtaje0wI4zxkxKq5mNA9mp8fiuR/vx3inVUpVBkt/TRXHcxn1cF8RnWO5N6JXRjVObJZr3/TWrM9o5ifdh/jp93HSLGWIDpR57ZBXD00kRmD4oh0wRpJfx8v3rpuMNPeWE1BeQ35ZZYNJXPvHIGft5fTzyeEcI/PNqSRWVQJQGSQL7ePde+mRkkGhceqrq3feu72cZ2adULSUvn7ePGn87tz7fAOvPbrAb7anG6bll97uIC1hwvoGh3M9SM7cPnAOEKaSUkfo8nMxiOF/LT7GD/vzuFYSdVJHxcV7Me0/u25fGAcfePCXP7/UGx4AP+9dhDXvb8Bk1mTnF7EM4t288KMfi49rxCiaVTWmJi9/JDt+L5JXdz+vijJoPBYi5Kz6rX2md6/aYt4tjbtwvz518x+3Da2E7OXH2ZBciYma1Z4MLeMpxfu5p9L9jKlZwyXDohlXLe2+Ps03WiX1ppDuWWsOpjPmkP5rE8pOOkaQLCs35naux2XDYxjdOfIs9oM4gwjO0fy5EU9eXbxHgC+2JhO37hwrh0uhdKFaOk+23CU/DJLCarYMP9m8XctyaDwSFrreqOCN46S1nNNpUt0MK9c2Z8Hp3Tl3VUpfLMlw5Z0VdeaWbIzmyU7s/H3MTCyUyQTukczNCmCbjHBTk26SqqM7MooZlt6EdvSjpOcXmR7Az6Z8EAfJveIYWrvmCZPVE/mltFJ7MgoYmFyFgB/W7SLHu1DGJTY8tZhCiEsqowm3lphL3V298QuzWIJiCSDwiNtOFLInuwSAPx9DFwrreeaXEJEIP+4tA+PTO3Ot9sy+WJjOnut/00Aqoxmftufx2/78wDLaFyf2DC6xATTKSqIDpFBRAX70jbEjxB/H/y8Dfh5GzBryxRvldFEflkN+WXV5JdVk1dazZH8cg7llnE4r4yckt9v/DhRfJsApvSM4fzeMQxLimjyEcDTUUrxrxn9OJBTxt7sEowmzd2fbuG7+8cQHSIbSoRoiT7fkEZ+meW9qV2oP1cOiXdzRBaSDAqP9OEae8PvGYPiCQ/0dWM0rVuIvw83jEzihpFJHMgpZWFyJj/sPEZKfv1NGhU1JjamFrIxtdBlsYT6ezOqcxRjukYxtmsUiRGBzXodaYCvF+9cP5hp/11NUYWRnJJq7v98G5/dNrxZJa5CiDOzjAoeth3fPaFzsxgVBEkGhQdKL6yo1+fxZinc22x0iwnhkak9eGRqD9IKKlh+IJd1hwvYllZ0yg0c58rXy0Dn6GAGJoYzMCGcgYlt6BQV1OKKYidEBPKfqwdy44cb0doy6v3Sz/t54sKe7g5NCHEW5m5Kt5Wqig7x46qhCW6OyE6SQeFxPlqbatvNOrZrFF1jQtwbkDipxMhA24ghQFZRJfuPlXI4r4zDeeVkF1fapn8rqk1U15qpMZlRypLo+XobiAjyJSrYj6hgXyKD/UiMCKRz22C6RAeT0CbAY0bPxnVry5+mdOMVa12yt1ekMDChDRf0aefmyIQQDVFda2L2cvuo4F3jO7t9XbIjSQaFRymvrmXuZnuR6VtGu7d2k2i42PAAYsMDmNgj+pSPMZt1ixvZc5Z7J3Zha9px2xrLR+Ztp3u7EDpGBbk5MiHEmXy1Kd02+xEV7NcsdhA78ozLZiGsvtmaQWmVpWdup6ggxndr6+aIhDO11kQQLL/7q1cNIL5NAACl1bXc/ekWKk9RHkcI0TwYTeYTRgU7NatRQZBkUHgQs7l+OZmbRie16uRBeJ7wQF9m/2GwrUzSvmOl/OXbnWit3RyZEOJUFiVnkVVsGRWMDPLlD8M7uDmi35NkUHiMFQfyOGLdoRri783MQc1jy74QztQ3Pox/TO9tO56/LZPPN6a5MSIhxKmYzZq3V9pHBW8enUSAb/MaFQRJBoUH+cChnMzVQxMI8pMlscIzXTU0gVmD7Rc7f1+0h+3pRW6MSAhxMsv25XIgpwyAIF8vrh+R5N6ATkGSQeERDuaUsupgPgAGhW2HqhCeSCnFs5f1oWf7UABqTGbu+Wwrx8tP3WFFCNH0HOsKXjMskbDA5tGb/USSDAqP8OHaVNvP5/WKISEi0H3BCNEE/H28eOu6QYT4W0bAM4sqeWBusq0ftBDCvTanFrL56HEAfLwUt45tvtUtmn0yqJQKUUpNV0o9q5T6QSmVr5TS1q8eDXi+QSl1h1JqnVKqSClVqpTappR6RCklbSk8QHGFkflbM2zHUk5GtBYdIoP495UDbMcrD+TxxrKDboxICFHHcVTwsgFxtA8LcGM0p9fsk0FgMrAQeAq4AIhs6BOVUj7Ad8DbwAggAPACBgAvAquVUsHODlg0rXlb0qkymgHo2T6UYR0j3ByREE3nvF4x3DOhs+349aUHWXso340RCSEO5JTy695c2/Gd4zu5MZozawnJIEAu8D3wd+COs3jec8BFQBVwExAIBAHTgEJgKJZEUbRQZrPm0/VHbcc3jOzQrHvNCuEKfzqvGyM6WS6CtIYH5iaTZ217JYRoeo6jguf3iqFLdPPuhNUSksHvtNYxWuuLtdbPAL805ElKqXbAA9bDx7TWH2mtTdpiMXCL9b5rlFL9nB+2aAqrDuWTWlABWMrJXDog1s0RCdH0vL0MvH71QCKDLCtf8kqreWhuMmZZPyhEk8surmRRcpbt+C6Hkfvmqtkng1rrcy2vPxPwA4qBd07yuguBA4ACrj3nAIVbfbIu1fbzrMEJBPpKORnROsWE+vPqVQOoGxhffSif/y0/5N6ghGiFPlp7lFrrhdiwpAgGJbZxc0Rn1uyTwUaYaP2+UmtddYrH/Gz9PqkJ4hFOll5YwdJ99jUZ149sflXdhWhK47q1rbd+8N+/HGBDSoEbIxKidSmvruXzDfalS7c14x3Ejjw5Gexl/b77NI/ZY/3eU8lCsxbnsw1p1HXhGts1io5RQe4NSIhm4KEp3RiaZBmJMGv445fbKCiT9YNCNIVvtmZQUlULQIfIQCb3jHFzRA3jyclge+v3rNM8pu6+YOvX7yilTvkl3KfKaGLuJnsLLikyLYSFt5eB/1wzkDbW4rY5JdU8PG+7rB8UwsXMZs0Hq+2dsG4Z3REvg3NzBVflJJ6cDNYNE1We5jEVDj9LiZkWZMmObI5XGAGICw9gUo9oN0ckRPPRPiyAV67sbztevj+Pd1aluDEiITzf0n25tg2Nof7eXOHQMrK58+Rk0Cm01qf8Eu7zsUM5mT+MSHT61ZcQLd2kHjHcMc5e2+yln/az5WihGyMSwrO953DBdc3wRIL8nL+h0VU5iScng+XW76cr+e3Ys6zMhbEIJ9qeXsT29CIAfL0MXDUkwc0RCdE8PTK1OwMTwwEwmTV//CKZkiqjm6MSwvPsyixmwxHLxZa3QXHTqCT3BnSWPDkZrFsPeLrCc3X3lWmtS10cj3CSj9fZRwUv6deeyGA/N0YjRPPl42XgjWsGEurQv/jpBbvcHJUQnud9h7WCF/Vt36xbz52MJyeDdTuFe5/mMXU7jve6OBbhJIXlNXy3w74nSMrJCHF68W0CeWGGva7+guQsFiZnujEiITxLTkkV3223fy61lHIyjjw5GfzN+n2sUsr/FI85z/p9aRPEI5xg3uZ0amotfYj7xoUxICHczREJ0fxd3K89MwfZF7M/9e0u0gsrTvMMIURDfbYhzVZkemhSG/rFt7zPJU9OBucD1UA4cNuJdyqlpgHdAQ180bShiXNhNms+32gvJ3O99CEWosGemd6LxAjLMunS6lr+9FUyJik3I0Sj1NSa+XyD/XPpplEtb1QQnJAMKqWClFKdlFLDlFLjlVJ9rH2BnUYpFVX3BTj2dQl3vE8pZft9tNbHgNethy8qpa5XSnlZX+8i4EPrfV9orXc4M17hGmsPF3DUoQ/xtH7Sh1iIhgrx9+HVqwbYdt5vSj3ObGlXJ0Sj/LArm3xrUfd2of6c37tlFJk+0Vkng0opL6XUpUqpN5VSO7D0/j0IrAOWAduBTKVUoVJqsVLqUaVUUiPjzHP42upw+7oT7ks84XlPAd9j2VH8MVCulCoHlgCRwCbgrkbGJprI5xvtG0dmDoonwNfLjdEI0fIM7tCGP07qajt+7deDJFt35gshzt6ctam2n/8wPBEfr5Y54drgqJVSHZRSLwGZWKZg7wb6WF9DneQrHLgIeAE4pJT6RSl1pXPDPz2ttRGYhiXhW49l2lgDycBjwBjZRdwy5JZW8fPuHNvxtcNPzPuFEA1x78TODO5gmWCpNWse/HIb5dW1bo5KiJZnR0YR29LsZc6uHtZyP5fOmAwqpdoqpf4D7AceBqKBHcD/gFuAwUAHIBTwBWKAHsAFwF+B77DU/JsMfKGU2mFdr9dgWmvVwK/UkzzXrLV+W2s9UmsdprUO1loP1Fq/qLWuOZs4hPvM25xRb4Fut5gQN0ckRMvk7WXg1SsHEGwtiJtaUMG/ftjn5qiEaHkcy5xd3K89bUNabpmzhpTHTsHS2u0I8BGWNXYHT/P4uinbA8DPANbdvBcD1wLTgQVKqUe01v9uROyilTCbNV84bByRUUEhGicxMpC/T+/Nw/O2A/DJ+qNc2Kcdo7pEuTkyIVqGgrJqFjmUk7mhhZc5a8g0cQZwE9BNa/2PMySCJ6W1rtJaf6O1ngn0xLJ5w/dsX0e0TisP5pFx3NJiOjzQhwv7tHdzREK0fDMGxTGlp32x+yNf76BMpouFaJC5DmXO+seHMTCxzRme0bw1JBnspbX+WGttcsYJtdaHtNa3Af/njNcTns9x2/7MQfH4+8jGESEaSynF8zP6EB7oA1i6kzz/vdTfF+JMak1mPnWYIr5hZJL7gnGSMyaDurHdj5v4dYVnySmpYum+XNvxNS14ga4QzU10iD9/n25v0vT5hjRWHcxzY0RCNH+/7s0lq7gKgMggXy7u1/Jnq85pD7RS6mnr1+gGPv5PSqmnz+VconWbuyndVhh3eMcIukQHuzkiITzL9P6xTHWojfbY1zsorTK6MSIhmreP16Xafr56WIJHzFada0GcZ4C/AUuVUnc04PGPWB8vRIOZzJovZeOIEC6llOK5y/rSxjpdnFVcxT+XyHSxECdzIKeUtYcLAPAyKP4wvGVvHKnT2OqIvsBsawHqlp+QALyRAAAgAElEQVQai2ZlxQH7UHxEkC8X9HFqYxshhFXbED/+cWkf2/GXm9JZcUCmi4U4keOo4Pm9YogND3BbLM7UmGQwB0vh6VosRZ1/VUpFOiUqIai/ceSKwfH4ecv1hhCuckm/9lzU137B9fg3srtYCEdl1bV8uzXTdnx9Cy8n46hRI4Na67eBKUABMA7YqJTqc/pnCXFm2cWVLJONI0I0GaUUz17ah4ggS9Wv7OIqXv5pv5ujEqL5WJScRXmNpbBKl+hgRnbynPGvRjfR01qvAoZg6UncEVirlLq8sa8rWrd5mzOw7hthVOdIOkYFuTcgIVqByGA//jatl+34o3Wp0rtYCEBrzWcb7OVkrh2WiFLKjRE5l1M6Kmut04DRwDwgGJinlJINI+KcmM2arzan245bcr9HIVqa6f1jGdetLQBawxPzd2I0md0clRDutSOjmN1ZJQD4eRuYOSjezRE5l1OSQQCtdaXW+irgKetNTyulvlZKBTrrHKJ1WJdSYOs4Ehbgw/m9Ys7wDCGEsyil+OdlffD3sXw87M0u4f3VR9wclRDu5biG/eJ+7Qmz7r73FE5LButorZ8HLgNKgcuBdYAkhKLB5m6yjwpePjDOI2o4CdGSJEQE8tCUbrbj1349QFpBhRsjEsJ9SqqM9foQe0o5GUdOTwYBtNaLgZHAYaAvEOKK8wjPU1xh5Mfdx2zHVw5JcGM0QrRet47pSK/2oQBUGc38ZcFOpHGUaI0WbMuk0mjZONKjXQiDEsPdHJHznWsyuBJYe7oHaK33AkOBn8/xHKIVWpCcaWv+3TcujF6xoW6OSIjWydvLwAsz+mKwrpFfdTCf73ceO/2ThPAwWms+W2+fIv7DcM/aOFLnnJJBrfUErfUVDXhcsdb6Aq21QWstc33itLTWfOkwRXzlUBkVFMKd+ieEc8PIJNvxc0v2UC61B0UrsjXtOPtzSgEI8PHi0oFxbo7INVwyTSzEudiVWcLebPturen9Y90ckRDiofO6ERVsrz345m+H3ByREE3nM4eNI9P7xxLq71kbR+pIMiiajbmb7X90F/VtT1iAZ/7RCdGShAX48PiFPW3H765KISWvzI0RCdE0iipqWLwj23b8hxGeW+bMuyEPUkotc8K5tNZ6shNeR3igKqOJhcn23VpXyRSxEM3GjIFxfL7hKFvTijCaNM98t4ePbh7qkWunhKjzzVb7GvY+caH0i/e8jSN1GpQMAhMADZzsL79ue9mZ3hVkG5o4pR92ZVNaZVmLlBQZyPCOEW6OSAhRx2BQ/OPSPkz772q0hpUH8vh5Tw5Te7c785OFaIFO7DjiieVkHDU0GfyYUydzVwF+wEdOiUi0So61BWcNSZARByGamT5xYfxheCKfWndW/uO7PYzr2pYAX9kbKDzPxiOFpOSVAxDs5+3xa9gblAxqrW861X1KqQuAaK31zc4KSrQuRwvKWZ9SCIBBwRWDPavNjxCe4s/nd2fJjmyOVxjJLKrkvVUp3D+5q7vDEsLpHAcoLh0QS5BfQ8fOWibZQCLczrEP8cTu0cSE+rsxGiHEqYQH+vLI1B6249krDpNbUuXGiIRwvuJKI9/vsm8cuXqo524cqSPJoHCrWpOZeZszbMdSW1CI5u2qoQn0aGdpKlVRY+Lfvxxwc0RCONei7VlUGS0bR3q2D6VPnOc3P5BkULjVigN55JZWAxAV7MekHtFujkgIcTpeBsWTF9lLzczdnG6rDyqEJ/jKYYr46qGtYw27JIPCrRzXZcwcFIePl/wvKUSdKqOJXZnF/Lgrm4XJmSzdm8PBnFLMZvcWZxjXrS3ju7UFQGv455K90rdYeITdWcXszCwGwNfbwGUDPLPjyIk8e0WkaNbySqtZti/XdjxriEwRC1FTa+bnPcf4eksG6w4XUG2tc+YoPNCHSd2juWpoAsM6Rrhl5OIvF/dk1cE8zBpWH8pn+f48JsrIvmjhHEcFL+jdjrDA1tH8QJJB4TYLkzOptY5wDOnQhi7RwW6OSAj30VqzIDmTV34+QMbxytM+tqjCyPxtmczflsmwpAgev6gHgxLbNFGkFt1iQrh6WCKfW9t1/fP7vYztGoW3jO6LFqrKaGKBQ/ODq1vRGnZJBoVbaK35eot948isIVJORjQfpVVGDuaWkXG8kszjleSUVFFSaaS40khFjQmT1mit8TIoAn29CfD1oiBpCl7GMuZuSiMpMohuMSG0CfJt0PnSCyt49OsdrEsp+N19HSID6dw2mAAfL4orjew7VkJ+WY3t/o2phcycvZYbRybx+IU98Pdpurp/D03pxqLkLMqqazmUW8aXm9K5boRnF+cVnuun3ccorjQCkBARwIhOkW6OqOk4ox1dRAMeA9KOTjjYnVXCvmOlAPj7GLiob3s3RyRaqyqjiS1Hj7P16HH2ZJewJ7uEowUVZ/9C7QYC8Ng3O203RQX70jU6hD5xoQzrGMmwpIjfTTutOZTPvZ9vpajCaLutTaAPN4xM4orB8SREBNZ7vNaanZnFfL4hja+3ZFBr1mgNc9amsvFIIW9dN5jEyPrPcZW2IX7cPaEzL/20H4D/LD3IzEHxUohatEiOa9ivHJyAweD5G0fqOKMdneNjTkdWFwubb7baRwWn9m5HiH/rWJch3E9rzZ7sElYeyGfNoXw2pRaedF2eM+SX1ZBfVsC6lALeXXUEpaB7TAgjOkUyuksUm1ILeX/1EUzW5RIGBbeN7cT9k7qc8m9CKUW/+HD6xYdzz4Qu/HXhLlYcyANgT3YJl/9vDe/fNJQBCU3TR/XWMR35eF0qOSXV5JZW8/G6VO4c37lJzi2Es6QVVLD2sGVk3qDgilY2W+WMdnRCnBWjycwih3UZ0nFEuJrWmn3HSlm8I4slO7JJPcPIn7dB0bltMB0iA4lrE0BsWADhgT6EBfgQ5OeNQSm8DIpak5mKGhPlNbX8/V+vUOsTzMSLL+dwXhmHcststcrsccC+Y6XsO1bKnLWp9e6LDvFj9nWDGdyh4Wv/EiMDmXPzUD7dkMaz3+2hxmSmoLyGa95Zzye3DmNIkut7fPv7eHH/pK48tWAXYClEfc3wRELlAk+0IPO22EcFx3drS/uwADdG0/Qa3Y5OiLO1fH8eBeWWNU/tQv0Z1TnKzREJT1VYXsPXW9KZuymdw9Y+oyfTuW0QozpH0Tc+jF7tQ+kaE4yf99lNdb56bCsAr1/9DwBMZk3m8Ur2HSth89HjbDhSyK7MYtsooKP+CeG8c/3gc+q+o5Ti+hEd6NU+hFs/2kxRhZFKo4mbP9zEZ7cPp1+860cIrxySwDsrU0grrKCowsh7K1P40/ndXX5eIZzBZNb1mh9c1Yo2jtSRDSSiyX3jsHHk8kFxeLWidRnC9bTWbDxSyOcb0/hh5zFqTL+fAg7282ZSj2jGdWvL6C6RLhkF8DIoEiMDSYwM5Pze7QAor65l89HjrNifx+6sYrSGkZ0juXtC50Zv/BjcIYKv7xrF1e+sI7+shtLqWm74YCNf3jGCHu1c20HB19vAQ+d15aG52wF4b/URbhiVRFSwn0vPK4QzrDyQxzFrW8WoYF8m9Yhxc0RN72w2kBRorWe5OB7h4Y6X17B0X47teOYgmSIWzlFrMrNkZzZvr0hhz0k6YgT6ejG5ZwyX9GvP+G5tm3TXbZ0gP2/GOxRsdrYu0cF8cutwrn5nPcWVRooqjNzw/kYW3z+GaBf3/J7eP463lqewP6eUihoTs5cf5q+X9HLpOYVwBseNIzMGxePr3frKIzX0N54AjHZhHKKVWLQ9C6PJMk3WPyFcaguKRqusMfHR2lQmvLycB75M/l0i2D8hnBdn9mPzU1N445qBTO3dzi2JYFPp2T6Uj28ZRrCf5Vo/t7SaOz/dQnWtyaXn9TIoHj6/m+34k/VHySo6fb1EIdyt8IQBiitb2caROq0v/RVu5biLWDaOiMaoMpr4YPURxr64jL8t2l2vULO/j4Frhyey+P4xLLx3NFcOTSDQt/WsiumfEM5b1w2mbgXGtrQinl6w2+Ut487rFWPbxVxTa+Y/Sw+69HxCNNai5EzbAMXAxHC6RIe4OSL3kGRQNJmDOaXsyLD2fPQyMK2f1BYUZ89oMvPFxjQmvrycfyzeU68Ac5tAHx6c0pW1j0/m+cv70icuzI2RuteYrlE8eVFP2/Hczel8au0W4ipKKR6dat848vWWDDKOn0PNRiGayDdbM20/t+ZlS63nUlm43dcOo4JTekUTHtiw7gxCAJjNmkXbs3jt1wO/Kw3TPsyfu8Z3ZtaQ+FY1Angmt47pyO6sEr7dZvnAe3bxHoYlRdC9netGP0Z1iWJYxwg2Himk1qyZvfww/7y8r8vOJ8S52n+slJ2Z1gEKbwPT+sW6OSL3kZFB0SRMZs2CbXIFJs7N+pQCpv13NQ/OTa6XCEYF+/L0Jb347c8TuHFUkiSCJ1BK8cKMvvRsb9lNXFNr5oEvt7l8/eADk7vafp63OYPsYlk7KJofx2VL5/WK+V13oNbkbN45w5RSHzTiXFprfWsjni9asNWH8skpqQYgKtiPcS7aTSk8y9GCcl74fh8/7j5W7/awAB/uHN+JmyQBPCN/Hy9ev3oA095YTXWtmX3HSnn5p/385WLX7fQd1TmSwR3asOXocWpMZt5ekcIz03u77HxCnK1ak9k2Yg5wRSsfoDibd1F/4MZzPI/C0sFEksFWyrG24GUDYvHxkkFpcWolVUb+u+wQc9ak1qsT6O9j4PaxnbhtbCfCAlrvVfzZ6hYTwhMX9uCZ7/YA8O6qI0zsHs2oLq4p+K6U4v5JXbjpw00AfL4xjXsmdHZ5eRshGmrVwXzySi0DFG1D/BjbtXU3PzibZNAIrHNVIMJzlVQZ+clhZGem7CIWp2A2a77ZmsG/fthn61JT5/KBcTwytTux4a2rTZSz3Dgqid/259n6GD82fwc/PzieAF/XlNkZ360t/ePD2J5RTE2tmbdXpkjdQdFsOK5hv3xgHN6tfIDibJLBQq31RJdFIjzWkh3ZVNdaRnd6tQ+1rV8SwtG+YyX8dcEuNqUer3f7oMRwnp7W21ayRJwbpRQvzerHlFdWUFJVS3phJa/+eqDejmNnn++Pk7ty60ebAfhsw1HuGt+ZtiHSlUS4V3GFkV92S/MDRx6fCiulDEqpm5VSvyql8pRSRqVUkVJqg1LqL0qp1llUqAk5ThHLqKA4UVl1Lc8t3sPF/1ldLxGMDfPnP9cM5Ju7R0ki6CTRIf485bBW8L1VKeyy7qZ0hUk9oukda7n4qzKaeW9VisvOJURDfbcjy7b8pG9cmEt317cUHp0MKqUCgV+AD4DJQBRQDoQCw4DngJ1KqU5uC9LDpeaXs/mo5QPe26C4dEDr3bov6tNas3hHFpNfWc57q49gMlsKv3obFHdP6MyvD49nev9YlJLe1c40a0g8IztFAmDW8Ng3O6g9Sf9mZ6gbHazzyfqjFFcYXXIuIRrKcRfxzEFxboyk+fDoZBD4KzAJy+aVJ4BwrXU4ls0w1wBFQAfgPbdF6OEWJNt3a03o3lYa1wsAUvLKuP79jdz3+TbbLnOAEZ0i+OGBsTx2QQ/ZJewiSimen9HX1n91d1YJH6076rLzndczhu4xlpGXihoTn25w3bmEOJPDeWVsSysCwMdLMX2AJIPg+cngtdbvH2qt/6W1LgbQWtdorb8EHrLeP1Ep1cYtEXowrevXFrx8oEwRt3aVNSZe+Xk/F7y2itWH8m23RwX78frVA/ji9hF0jZEpG1frGBVUrxbga78eIL+s+jTPOHcGg+LO8fbJlw/XHKHK6No6h0KciuOypUk9ookIkuYH4PnJYIz1+7ZT3L/F4edAF8fS6iSnF9kKBIf4eTO5Z7SbIxLutGxfDue9uoI3lh2yrdcxKLhpVBLL/jyeSwfEyZRwE7ptbEc6RgUBUFpVy8s/7XfZuab1j6V9mKWsTH5ZDfMdWoAJ0VRMZl2/tuDgBDdG07w0KBnUWhu01i1xsVeq9fvAU9w/2Po9R2st705OtjA5y/bzhX3b4e/jmhIWonnLLq7krk+2cMuczWQct3eiGJAQzqL7xvDM9N6E+kvNwKbm5+3F0w6lXuZuTmdHRpFLzuXjZeDWMR1tx++uSrGtERWiqaw9nE92cRUAkUG+TOguzQ/qnDEZdMXmCusO30Rnv+5JvGv9frNS6nGlVJj1/L5KqauAV7GsJ/xzE8TSqhhNZr7bbk8GL5N1Ga1Orcmye3TKKyvqdRAJC/DhhRl9mX/3KPrEhbkxQjGxRzSTe1hG7LWGvy3ajdlFSdrVwxIJ8besAz2SX84ve46d4RlCOJfjFPGlA+Kk+YGDhvxL7FNKfaSU6t7YkymlfJRSdwAHgZsa+3oN8BrwJpYOKC8ARUqpIqAS+BLYB0zXWn96qhdQSp3yS5za6oP5tqLB7UL9GW7dvShahy1HjzPtv2t4bsleymvs68OuGBzPsofHc82wRAwG+RtqDv56SS98rR+K29KKWLwz2yXnCfbz5voRHWzHb61IQWsZHRRNo7TKWO+idObgljlA4aqcpCHJ4EbgemC3Uuo3pdSdSqkGf7Iri4lKqbeBLGA20BbYfk4RnwWttQl4EHgYqLXeHIb99w6xxiKczHFdxqUDYvGSD/5Woaiihifm72Tm7LXszS6x3d41Opi5d4zg5Vn9iZQd5c1KUlQQN49Jsh2/+OM+qmtds8HjplFJtsQzOb3odwXGhXCV73dmU2W0rFXu0S6E3rEyK+HojLUbtNZjlFLTgeeB8cA44E2l1EEsGzB2APnAcaAGCAfaAB2BIVjW6wVhGZ0zYhmpe1Zrnef03+YESql2wEIsNQU/Av4NHAbaA1cATwMfKKW6aa2fONlryJXr2SurruVnhymgywa2zCsw0XBaa+ZvzeT57/fWayPn72PggcnduHVMR1spE9H83DuxC19tSud4hZGM45V8su4ot411fvnV6FB/ZgyK48tN6QC8veIwwzpGOP08Qpzom62OG0dabmWL0+UkjRkdbFAhL631IqXUd8AFwG3AJUB369c1p3lqXWQpWAo/f6i1ds0cxMl9jCURfF9rfZvD7YeAfymlMq2PeVQp9anWencTxuaxft59zHYF1j0mRNrPebhDuaX85dtdbDhSWO/2yT2ieWZ6bxIiZKN+cxfq78P9k7ryj8V7AHhj2SFmDU4gLND5G3tuH9eJuZvT0RqW7svlUG4pXaKlnJBwnYzjFWy0vj95GRSXyhr232nwpbq2+EFrPRPLyNosLGvy1mHZtVuGZeQvB9gNzAf+BAzTWnfRWj/flImgUqoXcJ718NWTPUZr/QlQgOXfYVoThebxHKeIZVTQc5VWGXnh+71c+PqqeolgbJg/b18/mPduHCKJYAty3YgOJFr/exVXGvnf8kMuOU/ntsFM7hFjO/5orRShFq7lWNliTJco6Y99EudU4l9rXQh8Y/1qrhy7rx85zeNSgEggyaXRtBK5pVWscSgmLO3nPI/ZrJm/LZP/+3EfeaX2QsVeBsWtYzrywOSuBPlJ95CWxtfbwKMXdOe+zy1lWT9cm8pNo5NoHxbg9HPdMjqJX/fmAJbWYH+e2p2wACkvJJzv980PZIDiZDx5EY9js83TlbGp295W6sJYWo3vtmdTV5liRKcIYsOd/0Ei3Gd7ehEzZq/lz/O210sEB3dow+L7x/DkRT0lEWzBLu7bnv7xloX1NbVm3vzNNaODIztH1mtRN29zukvOI8Se7BIO5pYBEODjxXm9Ys7wjNbJk5NBx93Kt5/sAUqpaUBdW4wNLo+oFXC8ApPagp4jr7SaR+Zt59I315Ccbi9MHBPqx2tXDeDru0bK2lAPoJTiT+fbq4jN3ZROxvEKl5znptFJtuM5a1OlCLVwCcfPpKm9Y+Ri9RQ8NhnUWqcAP1sPH1RKvaCUigZQSgUrpW4C5ljvTwUWNXWMnuZQbik7M4sB8PUycGHf9m6OSDRWldHE/5YfYtLLy5nnULDV18vAPRM6s+zhCVw2UNrIeZJxXaMY0sHSqt1o0ryx1DWjg5cNiCPcukEl43glS63TxkI4i8msWeTY/ECmiE/JY5NBq5uAvVh+z8eBHKVUCZYp4Q+BCCwbXmZorWtO9SKiYRZss//RTe4ZLWuAWjCTWfPV5nQmvrycF3/cT2l1re2+KT1j+PmhcTx6QQ+5yvZAltHBbrbjr7dmkJpf7vTzBPh6cfVQ+wqeD9ekOv0conVbn1JATollOUtUsC9jukS5OaLmy6OTQevu5cFYCk+vBAqBQKAE2Ao8C/TVWm9zW5AeQmvNgmTHQtNyBdYSaa35bV8uF72+ike/3mHr4wnQqW0Qc24eyns3DiEpKsiNUQpXG9U5ipHWrkEms+Y/Sw+65DzXj+xgK0i/LqWAfcdKzvAMIRrOsbLFJf1i8Zb2c6fk8f8yWutKrfXrWuvxWutIrbW31jpMaz1Ya/10UxS/bg22HD1OxvFKAEL9vZnYQxq7tDQbUgq45t313DxnE/tz7PupooL9eO6yPvz04DgmdI8+zSsIT/Kww+jgguRMDlkX4TtTXHgAU3vbF/TPkdFB4SRVRhM/7rI3P5BdxKfn8cmgaBqOV2AX94vFz9vLjdGIs7EhpYBr3lnPVe+sZ32KvV5goK8XD07pyopHJnDdiA7S1L2VGZIUwbhulos6s4bXfj3gkvPcPLqj7edvt2VyvFxW7IjG+3VvDmXW5S0do4LoFy/t505H3t1Fo9XUmlni0NxersBaBsckcF1Kge12L4PiuhGJrHhkIg9O6SbrAluxP51nHx1cvCObAznOr8A1pEMbesdadqJX15qZt0XKzIjGO7GyhWxyOz1JBkWjrTiQR1GFEbBM+9TtRBTNj8ms+Wn3MWa9tfakSeCswfEse3g8z13WV6r0CwYkhDOlp31pwOzlh51+DqUUN45Msh1/sTEds5SZEY1QWF7D8v32FWDS/ODMzumSXynlDZwPjAE6Wl9nH7BQa73ZeeGJluA7h6370wfEYjDIFVhzU1FTy9dbMvhg9RFSC+rXjfMyKGYMjOO+SV3oECkbQ0R9903qyq97cwFYtD2Lh6Z0IzHSuW0GL+nfnmeX7KG0qpYj+eWsSylgtOz8FOdoyc5saq0XFAMTw2XDWwOc6/xPGlC36rfuk18DTyqllgP3aq33NTI20QJU1NTyyx57fbDp/eUKrDlJL6zgi41pfLYhjeJKY737vA2KyyUJFGcwICGcMV2iWH0oH5NZ89bKwzx/eV+nniPQ15uZg+KZszYVgM82HJVkUJwzaT939s41GWwHbAPWAumAHzASmAJMBLYopW7XWn/ulChFs/Xr3lwqjSYAukQH06NdiJsjEjW1ZpbuzeHzjWmsPpSPPmHGLdTfm2uHd+CmUUm0C/N3T5CiRblnYmdWW3uOf705gz9O6ur0/3euHZ5oSwZ/3p1DbmkV0SHy/6c4O2kFFWw5ehywzHpcLM0PGuRck8ExWuu1J96olEoE3gQuBj5RSsVorV9tTICieVuU7DBF3D9WFum6idaa7RnFLErOYtH2TPLLfr8jMzEikFvHdOSKwfGyKUSclZGdIhmUGM7WtCJqTGbeXZXCXy/p5dRzdIsJYVhSBBtTC6k1a+ZtzuDeiV2ceg7h+RY61Lsd360tkcGy9rkhzukT4WSJoPX2NGCaUupO4HXgZaVUsNb62UbEKJqp4gojKw7k2o5lirhpaa3Zd6yU73dms2h7FkcLft9DVikY27Ut1w5L5LxeMVQZTeSXVbPvWCkFZdUUltdQUF5DUUUNlUYTlTVmqowmqowmlLIs7vdSCoMBgny9CQvwISzAh/BAH6KC/UiICCQxMpBQf+k248mUUtw3qQu3zLEsCf98Qxr3TuxCRJCvU8/zhxGJbEwttJ3jrvGdbUWphTgTrTXfOiSD0n6u4VwyPKC1flsptRP4FnhGKWXQWv/dFecS7vPT7mMYTZY5yH7xYbJItwlU15rYkFLI0r05/Lo3l8yiypM+LiLIl0GJ4cSHB3C80shbKw7z5Lc7KXRRDbfwQB86RAbROzaUfnFh9IkLo1tMCL7eUrDAU0zsHk3P9qHszS6h0mjiwzVHePj87k49xwV92tEm0IfjFUYyiypZeSCPiT2k0LlomJ2ZxaTkWVonBvl6cV7PmDM8Q9Rx2VyR1nqtUmo8sAx4WilVrLV+zVXnE03PsQH4tH4yKugKRpOZHRnFrE8pYH1KAZtTj9vWaJ7IoCy7uUzaUlqhbgdoUyiqMFJUUcT29CLqFgr7ehsYnNiGMV2jGNe1Lb1jQ2WneQumlOLeiZ2573NL9845a1O5fVwnp44K+3l7MWtIAu+sTAEsG0kkGRQNtWCb/TNpap92BPhK84OGcloyqJQKAOJO+IrFssGkHfASIMmgh8gtrWLtYcuCcqUspSFE4xhNZo4WlLMjo5hNqYVsSyvicF6ZbfT1TBpSms3Xy0DbED+ign2JDPYjIsiXyGBf2gT6EuTrhb+PFwG+XrYOMmatMZs1tWZNeXUtxZVGiiuNFFUaySmuIq2wgrTCCqprzb87V02tmXUpBaxLKeCln/bTJtCHid2jubhfe8Z2bSujhi3QhX3a06ntAVLyyimtqrVN5TrTNcMSbcngsn25ZBVVEhse4NRzCM9jMmsW77Ang5cNkCnis3GudQafxZLoOSZ+Z+r1Iu/8HuSHncdsycfQpAjah8mb9ZmYTGayiiupCE3CGBjJU9/uJP14JRnHK8grraa0qhZnlNoNC/ChQ2QgiRGBdIgMpENEEInW43ah/k4fnTObNXll1RzIKWVnZjG7MovZkVFs61Vd53iFkfnbMpm/LZNQf2+m9m7H9AGxjO4cJSOGLYSXQXHXuM48+s0OwNJL+JbRHZ2a2HeMCrKVsjFrmLspnYccOqEIcTIbjhSQW1oNQFSwL6M6R7o5opblXEcG/4KlruCp3sFrgOhFXugAACAASURBVGwgE8iyfs88xWNFC1RvirgVbRw5kFPK6oP5bM8oIrWgnLIqy2aLGpMZs1lbRtK0ZURNa02tSWM0W0bXbIler1kAfLoh7ZxiUArah/rbErwOkUH1Er+wwKbdzGEwKGJC/YkJ9Wds17a223NKqlh9MJ9VB/NYfSi/3g7nkqpa5m3JYN6WDBIiArh6aCKzhsRLKZEW4NKBsbz0837ySqs5VlLFd9uzmDk43qnnuHZ4or2UzZYMHpjcVS4YxGl9t93eEvWivu3xll7qZ+Vck8HtQAb2RM8x4cvSWuc7JzzRHGUcr1/H6aI+7dwckesUVdQwe/khftufx5H88gZP2TaWr5eB6FA/OrcNonPbEBIiAqyjfUHEtwnA36f5r4WJCfVn5uB4Zg6Ox2zW7MoqZsnObBZvz6638SW9sJKXftrPq78c4LxeMdwypiNDOrSRMkXNlJ+3FzeNSuKln/YD8O6qFGYMcm7v18k9o+ttJFl7uIAxXaUItTi5mlozP+yyJ4OtaYDCWc61tMxAZwciWo7FO+x/dGO6RHlcHSez2cyctUf5aF3qScu1OOEEKLOR6IgwIoJ8iAsPoFNUED3bh9I1JoSENoGEBnh7VDJkMCj6xYfTLz6cxy/oQXJ6EQuTs/h2W6atM0qtWfPDrmP8sOsYAxLCuWNcJ6b2bielRZqh64Z34M3fDlFRY2LfsVJWHsxnfLe2Z35iA/l5e3HZwDg+XJMKwFeb0yUZFKe05lA+RRWW95HYMH8GJ7Zxc0Qtj1SeFWftxELTnqKqppbnf9jHV5vTqTL+fkNEnQAfL6KCfUmICCQy2JdQfx+C/SybLry9DPh4KXy8DHgbFCH+PrQL87d8hfoR5OfDhAkTAFi+fHnT/GLNjFKKgYltGJjYhscv7MH3O7P5YmMam1KP2x6TnF7EPZ9tJSEigDvHdWbWkHjbphbhfmGBPlw1NMGWrL2z8rBTk0GAWYPtr//j7mMUVxibfAmEaBkcly1d0j9WlhScA0kGxVk5lFvGnuwSwFI65PzeLb+Ok9ls5v9+3Mf7q1Ntzc0deRsUAxLCmdonhisHJxAW6NxCu62Zv48XMwbFM2NQPPuPlfLB6iN8uy2TGpMlGU8vrOSpBbt487dD3D2hM1cOSWgRU+StwS2jO/LxuqOYzJo1hwrYlVlMn7gz7SNsuF6xofSNC2NnZjE1tWYWbc/k+pFJTnt94RmqjCZ+3n3MduxJAxRNSVZYirPieAU2qXs0IS2888SCbZn0+/svvL3yyO8SwS7Rwbx4RT8OPHcBX989itvHdpZE0IW6twvh/67ox+rHJ3L/pC6EO4wCZRdX8fTC3Yx/6TfmrDlC1SlqLYqmkxARyEUOfV/fXZXi9HNcOcS+MeWrzRlOf33R8i3bl0t5jeX9oGOUpfC9OHuSDIoG01qz2CEZnD6g5V6BFVXUcOmbq3lwbjJl1bX17hveMYJfHxrHr38az5VDEjAY5M+kKUWH+PPw+d1Z+/gknrq4J1EOa1JzSqp55rs9jHtRksLm4I6xnWw/L96RfcqOOOdqev84W9manZnF7Mkqcerri5bvuxMqW3jSWuumJJ9yosF2Z5WQkm9v9TOphXYG+HxDGkOe+5Xt6cX1bu/dPpTfHh7P3DtH0iUmxE3RiTqBvt7cNrYTqx6dyF8v6UXbEHtSmFtqSQonvrycT9cfpeYkRa+F6/WND2NkJ0s9N5NZ8+HqI059/bBAHy7oba9WMG9LulNfX7RspVVGlu6zd1qaLs0Pzpkkg6LBHKeIz+/drsWt3aqtNXPThxt58tud9aaEQ/29+eTWYSx5YCwd2wa7MUJxMgG+Xtw6piOrHp3IM9N6ERNqTwqzi6t4asEuJr68nLmb0jCaJClsareP62j7ee7mdMpPGGlvrCuHJNh+XrAtk+paGQ0WFr/sybFdCPZsH0qXaLmIP1eSDIoGMZt1veH4lrZI90heGcNeWMry/Xn1br98YBzJT59Xr1iyaJ78fby4aXRHVjwykacv6VVv+jizqJLHvtnJlH+v4JstGdRKUthkJnSLplNUEAClVbV8s9W5a/tGdY4kztqO7niFkaVN2HNbNG/1mx/IqGBjSDIoGmRL2nGyi6sACA/0aVE1v5buzeG8V1dSWG7vgBHo68XXd43k1asGyJrAFsbfx4tbrCOFT17Ug4gg+6aeowUVPDxvO+e/tpKFyZmYGtKwWTSKwaC4cVSS7XjOmlTMTvx3NxhUvQ4n8zbLVLGAwvIaVh+097eY1q9lDVA0N/IpKBpkiUOh6Qv7tMOnhbT6+e+yg9z60eZ608L948PY/NQUhiRFuDEy0VgBvl7cMa4zqx6dyCNTuxMWYN99nJJXzgNfJnPh6yv5fme2U5MT8XszB8cT4mepVJaSX86Kg3lneMbZmeWQDK48mE+etQetaL1+2JVte18fmBhOQkSgmyNq2VrGJ7pwK7NZ8/1OezJ4SQu5Anvk6+28/POBerfdOa4jC+8bQ6CvlNj0FEF+3tw7sQurH5vIQ1O6EeJv/297IKeMez7bysVvrObn3cfQWpJCVwj28+bKofa1fXXFop0lISKQoUmWrhIms2bxjqwzPEN4upa8bKk5kmRQnNHmo8fJtV6JRwb5Mrxj8x9Ru3XOJuY51CUzKHj96gE8cVEvN0YlXCnE34cHpnRl9aOTuH9SF4J87Ruc9maXcMcnW5j+3zX8sDNbpo9d4MaRSdRV9Vh5II9DuWVOff3LBsbZfv52W6ZTX1u0LMeKq9hwpBAApeDivrJesLEkGRRntMThKvyCPu3wbsZTxGazmZmz19QrN+DvbWDx/WO4dEDcaZ4pPEVYoA8Pn9+dVY9N4q7xnQlw2PW+M7OYuz/bysSXl/PxulQqapy787U1S4wMZEpPe0eiOWudW2bmkr6x+Frfe3ZkFDs92RQtx5Kd2dQN8o/oGEl0qL97A/IAzfdTXTQLJrPm+132Vj8X92u+V2Bms5nL/7eWLUeLbLeF+Huz9M/j6RXrvDZZomWICPLl8Qt7sOqxidw2piN+3va3u7TCCp5euJtR/1rGC9/v5Yi1fqZonJtHJ9l+/mZLJsUVRqe9dligT73apgtkdLDVWuQhzQ+aE0kGxWltSi20LdaOCvZleMdIN0d0cmazmRmz17E9w15IOirYl1WPTiQuXBYWt2ZRwX48dUkvVj02kfsmdqm30aSowsjbK1OY+PJyrnlnPYu2Z0kdu0YY2SmSHu0std4qjSbmbk5z6uufOFUsG4Nan7SCCranWy74vQ2qXlFyce4kGRSnVX8XcXu8DM2z1c9V76wnOd0+Ihgd4sfyP08kXHoJC6voEH/+PLU7656YxN+n9yYhIqDe/etSCvjjF9sY8uyvPDQ3mV/25Ei7u/9n77zDojrT93+/U2CAoTcRUQQpYkOxx67RxDRjEtPbpm/6JtndfONukk3bJL9N2ZLd9F6MaZpoNPbesCtFiggivbdhGOb8/jjjKUgZmDP9+VyXV+aVc868BpjznKfcdz9hjMmyg5/uPqOo5uPc1EghmC+tb0PmmTrFrk24Bz9L2pZmJUciNIA+45WAgkGiRzrNHH6VlIgXu2iT7v2fZ+JAkXhTiND7YvMTc6DX0cQwcSH+PhrcPj0eW5+ci4/umIgFI6MhfcZpajfhx8OluOezTGS8sAEPfnkIX+0rxpmaFsWmkc1mDhWNBpQ3GDxuwvmq9FiE+osB28bsCsWu7atRy1pVfjysrMA14fr8IklQkNC0ctDdkuiR/adrUd18vkTsi8kuOEX83OqTWHdSvNmEB/hg65MUCBJ9o1YxzEuNxrzUaJQ1tGFl5lmsPFiCkto24ZgWYyfWHC/DGou0UlSgL8YOCcbo2GCMiNIjPjwAcaH+CPLr/ueN4zi0GDtRXNOKwupm5FU041BxHY6U1KPJwA+vJEQE4O0bxmPMEM/oa9Vp1bhpylD8Z0sBAOCjXUW4ZLRyN+2l42Px1T6+/PzLsTI8e8Uot7PGJAZGYVUzsssaAQA+GpVsYImwDbpjEj2y5riYjl88ZpDLlYjf316AT3YXCesAXzU2PD6bAkGi38QE++GR+Ul4eN4InDzXiLXHy7D2eBmKalplx1U2tWNjdiU2drFE06gYzBm/h6qzA7Nf3wKO43vmGlo7YOyjTFpY3YLr3t2N168dhys8RC/t1qnxeHdbIUxmDvtP1yLrXCPSBgcpcu2MYaGIC/NDSW0bmgwmbMmpxKUuWrUglEWqdzs7ORKBOm0vRxP9gcrERLeYOs1YJ50idrEP21VHSvHS2hxhrVUzrH1kJsL01D9CDBzGGEbHBuOPl6Riy5NzsP6xWVh+2UjMS42S6RZ2xWTmYNYGwKQLwZmaVhTXtqKqqb3XQDBI8tBi6DDj4a8P4x+/5XrEUMSgYB0uGS029n++94xi12aM4WqJTNQPNFXsNaw57rr3JHeHUihEt/AlYt7LNyrQ16Ws2w4X1+HxFUeEtYoB39w7DcPCA5y4K8LTYIwhZVAgUgYF4u6ZCeg0cyisasbx0gacPNeIMzUtOFPTitL6NrQaex400WlVGBzsh+ERARgeEYDRscHIGBaKIaF+OF3dgrs/zUShRdrmX5vzkVvehLduSHd7l5zbpsUL/V0/HS7F04tTEaRQJmfJ+Fj8c3M+AGBrbiUaWjsQ7E9ZIk+ma4l4/sioPs4g+oN7f9oQduMXSTp+8RjXmSKubTbixvf3Qpo8+c/NE5AxLNR5myK8ArWKISk6EEnRgVg6Qf61dlMn5l96JTi1Fl9+8SUAwFerQqi/T6/9bAmRevz44EV4+OvD2H6K9/P9LasCN763Fx/eMQkRel+7/XvszaT4UKQOCkROeRPaOjrx/cGzuPOi4YpcOyFSj7FDgnHsbAM6OjmszyrHsolxfZ9IuC1UIrYvVCYmLuCCErGLCE2bzWZc8e+dMHSIpbdnLkvFpQo2pxPEQPDVqKHpaIHWUI/4iADERwQgJtjPqsGGYD8tPrp9Iu6eIQZKR8824Jr/7kaRG4thM8Zwy9RhwvrzvWcUnZy+QuKRLvWpJTwT6RTx5S5yT/IkKBgkLmDf6VrUtvAl4uggX2QMdY2s2+8+zURpvTjpefnYGNwzM9GJOyIIZdCoVVh+eRpeWDJakLk5U9OKa/67W6af6W4sGR8LvS9fgCqsasGeghrFri19SN1dUCMoHxCeR0FVM3LKmwCcLxHTFLHSUDBIXMAvXYSmVS5QIn5r4ylsza0S1omRevzzhnQn7ogglOfWqcPw31syBOu8mhYjbnxvL7bmVvZxpmui99XgmgnisMdne5QbJBkc4odJ8fyDaldNVMKzWCu5J81JjhQeMAjloGCQkMGXiF0rHb8jrwpvbcwT1npfDVY9eBFUKvrxJTyPRaMG4at7piDEMhDR1tGJez7LlP1euhPSUvGG7AqUNbT1cnT/kErxUKnYc1kj6Rd0lbYlT4PupoSMPYU1qLOYyw8K0mGCk0vEDa1G3P1pprBWM+CH308nLUHCo8kYFobvH5iO2BDeMq+jk8ODXx12S8eNpOhATEvgPc07zRy+3qecX/Glo2OEsvqBolpFA03CNaASsWOgYJCQIfUiXjzG+SXia/+3B+0mcWDkH8vSkRwd6MQdEYRjSIzUY+X905AQwUsmdZo5/OHbo/hyn3KlVkdx6zQxO/j1gRIYTcr4FUcG+mJ6YgQAgOPkn1+EZ0AlYsdAwSAh0NFpxrqTrjNF/Nzqk8irbBbWV4yLwZLxsb2cQRCexeAQP6y4bxpSB/EPQBwHPPPjCXywo9DJO+sfF6dFIzqIl8mpamrH+pPK9fdJ/Wl/pmDQ46ASsWOgYJAQ2F1Qg3pLiTgmWIfxcSFO28vOvCqZ1VxMsA5vX08DI4T3ERnoi2/unYpxEu/iF9dk46Odp524q/6hVatw4+ShwlpJR5JFowZBq+YrGEdL6lHcxUKQcF/yK8USsS+ViO0KBYOEwFoXKRE3tBpx92din6BGxfDdA9NoYITwWkL8ffDF3VOE6VkA+NsvWfhsT5HT9tRfbpw8FBrLZ8r+07XItdzkbSXE3wezkiKF9c/HaJDEU5AKTc9JoRKxPfGauytjLIUx9i/GWC5jrIUx1sAYy2aMfcQYm+3s/TkbVyoR3/D+Xpmw9OvXjkNsiL/T9kMQrkCgTouP75yMiRK3nb+uOokvFMyy2ZPoIB0WjZL6FRcpdm2aKvZM1spKxIN7OZKwFa8IBhljjwA4BuAhAMkAzAB8AKQCuBPArc7bnWuwK78aDW18iTg2xM9pJeJ3tuYju0zMGCweMwhXT6A+QYIAeFmlj++chAlDxd/P5T+dwDf7lZvQtSdSmZmfDp9DS7tJkesuSIsWtBlzypuQL+k1JtyTC0rEqeRFbE88PhhkjN0H4G3wPsyvAhjGcVwgx3F+AGIA3AZgtxO36BLIp4gHgTHHl4hLalrx/9bnCuuIAB/8+8bxDt8HQbgygTotPvndZIyTPLA9/eNxrHaDjNjUhDAkRvLT0c3tJsWyeHpfDeamiMGCu2oyEiLSrODclCgEUInYrnh0MMgYiwfwhmV5P8dxf+Y4TniE5jiunOO4zzmO+8gZ+3MVOjrN+C2rQlgvHuOcEvEN7++F2WJdygB8cfcU6hMkiG4I0mnx2e8mY6xlqITjgD+sOIItOa7tVMIYkw2SfKVgRvPSMWIJeu1xciNxd2QJCpoitjuefqd9FIA/gH0cx73v7M24KnsLa4QS8eBgHdKdUCJ+dvVJme/w/bMTkRoT5PB9EIS7EOynxad3TkZSlB4AYDJzeODLg9h/utbJO+udayYMgY+lpHvsbANOlDYoct15qVHCdbPKGnGmpkWR6xKOJ7+yCbkVVCJ2JJ4eDN5k+e/XTt2FiyP19Fw02vEl4mNn6/GpREYmPtwff7o01aF7IAh3JDTAB5/fNQVDQnmnEkOHGXd9ckCxAMsehAb4YPFoMYunVHYwUKeVTRWTV7H7suaY+L2jErFj8NhgkDGWCOD848RhxthUxtjPjLEaxlgbYyyHMfY6Y8yrHzk6zRx+k0wRXzrasel4s9mM2z/aL6zVKoav75nq0D0QhDszKFiHL+6aggg9L+rc1G7C7R/td+nM2E1TxEGSVYdL0azQIMmlkiDz1+PUN+iurCWhaYfjscEggCTJ6zkAdgK4HIAWAAcgBcCTAI4wxkb1dBHGWI9/PIHMolpUNxsBABF6X2QMc6wX8Z9/OC54IQPAXy9PQ4zFj5UgCOuIjwjA53dNRpDFs7umxYg7Pj6A2hajk3fWPZPiQzHCUt5uMXZi9RFlBkkWjIwWBajPNuBsHQlQuxtdS8TzqEQsw14xiScHg9LGt2cBnAIwleO4IAB6AIsBVIKfKP6eMeaVeWhZiXhUNNQOFJrOKmvAt5lnhfWowUG4fXq8w96fIDyJkTFB+OiOSULf3OnqFtzzWSYMHZ1O3tmFdB0k+VqhUnGwvxYXjYgQ1uuoVOx2SEvE81KpROwoPDkYlP7bOABXcxy3DwA4jjNzHPcrgN9Zvp4CYGl3F+E4rsc/7o7ZzMk+LB09RXznRweE1xoVw6e/m+TQ9ycIT2NifBjevj4d55MEB8/U4fEVR2A2u97n1TUTYoXA9XhpA46fVabPcbGk1WUtlYrdDun3zFnKFq6MvWISTw4Gpaqj6ziOy+16AMdxa8BnDAFgvkN25UIcOVuP8kYDACDUX4spw8Mc9t4vrslCRVO7sP7TJSmI0Osc9v4E4alcOiYGyy9LE9a/nijHy2uznbij7gnx98Hlkpv9V/uVcVK5OE2scBwqrkdZQ1sfZxCuQkFVM5WInYQnB4PSJpQLAsFuvhZnx724JNKs4MVp0dCoHfPjcKamBR/uPC2sEyMDcM+sRIe8N0F4A3fNGI47L4oX1h/sPI2Pd53u+QQnceMUsVS86sg5NBk6ejnaOkIDfDA9MVxYU6nYfZB+r2YnR1KJ2IF4cjCYBd52zlpcr45iRziOw68SlX5HThHf+uE+nM9oqxjwye8mO+y9CcJbWH5ZGhaNihbWf/sly+UCo4nDQgWdxFZjJ1YpNEgi/Tz7lQSo3Qbpz6dURJywPx4bDHIc1wpgj2WZ0suh579WZNcNuRgnzzWipJYvnwT6ajB9RHgfZyjD+9sLUFwrlm0emDMCcaH+DnlvgvAm1CqGt28Yj/EWH2OOAx5fcQQnz7mOBiFjDDdJsoNf7StWpB974ahonJ+FO3CmFpWWdhjCdSmpbcVxiz6mVs0wLzW6jzMIJfHYYNDCZ5b/XsIYuyAgZIxdBiDZslzrsF25ANKsIG/yrrb7ezYZOvCaxHs4JliHpxb1FqcTBGELOq0aH9w2EcPC+Qeuto5O3PvZQVQ3t/dxpuNYOn4IfCXOIccUGCSJ0PtisqUHmuOADdkVfZxBOJv1Er3bi0ZEINhP68TdeB+eHgx+BL5crAbwA2NsMgAwxlSMsUsAfGg5bi+8KBjkS8TiL94lox2Tjr//84Po6BSf+t+/baJD3pcgvJlwvS8+uG0i9Jb+q9L6NjzwxUEYTf3porEfwf5ambDwV/uUkZm5ZJT4ufbbSQoGXR3pPelSB92TCBGPDgY5jjMBuAJACYA0APsYY40AmgD8CiAafLB4LecJWjFWklfZjMIq3p3A30eN2cmRfZxhO7vzq7GroEZYXzp6EEbHBtv9fQmCAJKiA/GvG8cLkjMHiurwl59OuIxE1s2SUvHqo+fQqMAgyYI0scy4u6BakeEUwj5UNhpw8EwdAL6P/OI0CgYdjUcHgwDAcVwhgDEAXgIf+GnAD4scAvA0gMkcx5U6b4eOR9pQPTclCjqtfUvEZrMZD351SFj7adV4+4bxdn1PgiDkzE2Nwp8vET2/V2SW4BOJJ7gzmTA0FCnRgQD4UrYSgyRDQv0xanAQAKCjk8O2U1U2X5OwD9IS8ZTh4QgL8HHibrwTjw8GAYDjuAaO45ZzHDeK4zh/juP0HMdlcBz3d47jXNfA005I+wUdUSJ+aW2OzHLuxSWjBLFZgiAcx72zErB0fKywfuGXLOzIc36Q1HWQ5GuFSsUL06hU7A78SlPETofuyF7G6eoW5JTzop4+GhXm2lnUs6LRINM3S4kOxDUZXifpSBAuAWMMLy8dg/Q4fsLYzAEPfnkIRdXOfyZeMj5WNkhyotT2QZKFEmmdLTmVLtMnSYjUthix73StsF40ioJBZ0DBoJchzQrOTo4Umsrtxb2fZeK8ExZjwAd30NAIQTgTnVaN927NwKAg3vGn0WDC/V8cRJvRuR7GwX5aXCZxJPnmgO3ZwdRBgRgS6gcAaGo3YW9hTR9nEI5mQ1Y5Oi03iYxhoYgOIicqZ0DBoJexzoETW9tyK3FUIhNx0+ShpClIEC5AVJAO792WAR+L61BOeROe+fG40wdKlk0SqwarDp+zOUBljMlLxVkkQO1q0BSxa0DBoBdxtq5V0PDSqhnmj7SvqOfj3x4VXgf4qvHCVaPs+n4EQVjP2CEheF7yO/nD4VJ8oVCv3kCZMjwM8RZNxKZ2E9YeL+vjjL6Rloo3ZlXCbHaNCWoCaDR0YFd+tbCmErHzoGDQi5BmBacn2lfU880NuahtMQrrl5aMgUpFP24E4UrcMCkO12UMEdZ/+/kkDhfXOW0/jDFcP0kcJFlxoMTma04cFopQf/6zrrzRILhcEM5nc3aloD07OjYIcWFUOXIWdHf2IhxVIm4ydOA/WwqEdWJkAJZIJhgJgnANGGN4YclomQTL7788hBonOpRckxELtcVLbn9RLQqrmm26nkatklmbUanYdZD2sEv9pAnHQ8Ggl1DRaMBByxO/WsWw0I7p+Ie/PgyTpBTz7i0ZdnsvgiBsQ6dV4783ZyBIxw+TlTUY8Og3R4SmfkcTFajDfInKwYpM27OD0lIxScy4Bq1Gk0z70VFOWET3UDDoJaw/WY7zveFThofZTdQz61wDtuaKv+CLRkVjhEVMliAI12RouD/euiFdWO/Mr8YbG3J7OcO+XC8ZJPn+4Fl0dNomCTMrKRI6LX+7y6tsxmkXkNLxdrbmVsHQwX9fk6P1SIzUO3lH3g0Fg16C1HXEniViqdOIj1qFN69P7+VogiBchXmp0Xhk3ghh/Z8tBU5z7ZidHInoIF8AQHWzEZuyK226np+PGjOTRNvNDVQqdjrStqVLqETsdCgY9AJ4UU9eX4sx+01srT1WhtPVrcL68YuT4O9jXx1DgiCU49EFyZiZFCGs/7DiCCobDQ7fh0atwnUScfoVCmgOLkyjUrGr0G7qxOYcMcAnSRnnQ8GgF7Axq0IQfs4YGoooO4h6ms1m/N9Px4V1WIAPHpgzopczCIJwNdQqhjevT0dkIJ+Vq2kx4rEVzukfXDZRDAa3napCWUObTdebPzIalrkUHCyuQ1WT84ZkvJ2dedVobjcBAOLD/ZE6iFqJnA0Fg16A1ATcXlnBf28pQL3Ef/jvS8fY5X0IgrAvEXpfvHV9OpglcNpdUIN3tuQ7fB9Dw/1x0YhwALxt3neZZ226XliADybGhwEAOA7YlE3ZQWchFZpeNHoQ2PkfNsJpUDDo4TS3m7DDzqKeBqMJ/94s3ixGROrtOq1MEIR9uWhEBB6eK2b239x4Cvsl/rGOQpodXJFZYrNgtKxUnEXBoDPo6DRjg+T/PUnKuAYUDHo423KrBHP21EGBGBquvKjn0z+egFEy7ffPm2hohCDcnUfmJ2GyJZNm5oBHvj4sE5J3BItGDRLE8c/WtWF3gW3ewlJrup351Wg1mmy6HtF/9hXWoqGNryINDtZh3JBgJ++IACgY9HjsXSKubDTgpyOlwvqixHCkxdAvN0G4Oxq1Cm/fmC5z73hq5VGH+hfrtGpcLRGs/8bGQZKh4f5IiuIlTIwmM3bn2xZcEv1HKjRNJWLXgYJBD8ZoMmOLZGLLHsHgQ18fFvQLVQz4140TFH8PgiCcQ0yw8IHxAAAAIABJREFUH/6xbJyw3pRTiY92FTl0D1LNwd9OVqDOxuzkvJGioPWmHNska4j+0WnmsP4klYhdEQoGPZjdBdVoskxsxYX5YWSMshNbx87Wy/qIrpsYhzC9fcSsCYJwDvNSo3H3jOHC+tV1Ocgpb3TY+4+MCcK4uBAAgLHTjB8Pl/ZxRu/Ml1jTbcmpdGim09s5VFyHaovVYYTeFxnDQp28I+I8FAx6MNInsEVpyqfjH/n6sPBap1HhhatGK3p9giBcgz9ekir4FxtNZjz2zREYOjod9v43TJJqDpbYFMBNGBoi9CGWNxqQVea4wNbbWS+ZIr44LVrwoCacDwWDHkqnmZNNbC1SWNRz/YlyFNVIBaaT4aOhHyeC8ER8NCq8fUM6fC2/4znlTXh9vePs6i4fGwM/rRoAkFvRhCMl9QO+lkatwpwU0Y1ks43uJoR1cBwnm+BeJPGLJvqP2WxGVlkDfjxUiv9tLcCLa7Jsuh7ZQ3goh2XpeB9MGKpsOn75qhPC67AAH9w3O1HR6xME4VqMiArE8stG4i+rTgIAPtx5GnNSImU2b/YiUKfF5WNjsPIgrzW44kAJxtvwmTYvNQqrjpwDwPcNPjw/SZF9Ej2TU96E4lo+gRDoq8H0xIg+ziDOU91swHeZpdhTWIO8iiZUNbejo1PZ9gYKBj0U6RSx0un4r/YVy9T7n7syTbFrEwThutwydRi25FYJVmJPrjyKdY/OQmiA/XuFb5gcJwSDq4+ew18uT0OA78BuYbOTI6FWMXSaORw9W4+qpnbBdYWwD1ILwDmpUVRJ6gWz2Yy1x8vxxb4zOFHaKLi12BP6bnggHCef2FJSANpsNuPvv2YL69gQP1w5LraXMwiC8BQYY3j1mrEItwR/FY3tePqH4w4ZwpgwNBQjLLIwrcZOrDle1scZPRPi74MMS2aR44CtuVQqtjdymTMqEXfH1txKLPnPLiQtX4eHvj6MvYW1fQaCvhoVgnQaRNn4MEOZQQ9Emo7X+2owPTFcsWv/d1shGg3iD+er15DtHEF4E5GBvnjt2rG469NMAMC6k+VYefCszC3EHjDGsGziELy8NgcAb09ny3vOGxmF/UW8GsLmnEpcZ+f9ezMlta3CoI6PWoXZyfZvLXAXmg0mvLQmC6uPnUNLe89DWT4aFYaG+WPU4CBcNCIC0xLDERcqN5Fgywe+DwoGPRDpE9iclEj4atSKXNdkMstt56L0mOGAfiGCIFyL+SOjccvUofhiLy8C/fzqk5g6PNwuDkdSloyPxavrctFp5rC/qBaFVc1IiNQP6FrzU6Pw91/5wHJHXjWMJjOVLu2EdHDkohHhCNRpnbgb16CkrhX/98Nx7MqvRk8ui7Ehfrg4LRo3TRmK5GhlpeG6Qj/5HohMUkbBEvEr63LQJpGTeEMiRksQhHfxzOI0JEQGAABajJ146rujNnsH90VUoA5zU0TR6O8sPYQDYUSUHnFhfgB4D/cDRY73XvYW7O2E5U5UNhpw64f7MPPVLdiRd2EgGBXoiwfnjsDx5xZi15/n4bkrR9k9EAQoGPQ4SmpbkS1Jx0slFGzBYDThsz1FwnpcXDDGDglR5NoEQbgffj5qvLksXRhO23e6Fp/sLrL7+y6bOER4/f2hs+gcYADKGJMJUG8iiRm7UNPcjkxLoM0Yn1X2Rhpajbj3s0xMeWUTduRVy76mYsDMpAj89vgs7H9mAZ5alOLw7CkFgx6G9AlMyXT8Mz+dkI2y//OG8YpclyAI92VcXAh+P0eUlXp1XQ4Kqprt+p5zU6MQoRcHWLbnVQ34WvNSpdZ0FeRGYgc2ZVcK2a+Jw0K9cmr77U2nkPHiRvyWVQHpj5hGxXD9pDgce3YRPr9rikMygD1BwaCH8ZsdSsQNrUaZBdSMEeEYFh6gyLUJgnBvHp6XhJExvDtJu8mMJ1cehanTbLf306pVuHq8qGCwMrNkwNeakhAGfx++p/pMTSsKq1ts3h8hR5qgWJjmXSXiw8V1mPLyRry5IQ8mSQZbzYCl42Nx7LmFePWasdDrnD++QcGgB1Hd3I4DZ/h0vIoBC9KUScc/ufKo8GTHGPDGsnRFrksQhPvjo1HhH9eNg1bNl4sPF9fjvR2Fdn1P6eTvhqwK1LYYB3QdX40aM5NE8WNyI1GWlnYTduSLJVFv6Rc0mcy459MDuPqd3ahobJd9bVFaNI48uxBvXJ8Ofx/nB4HnoWDQg9goSUFPHBaGCL3t6fjqZgM2Sj4gLxk1CFFBOpuvSxCE55A2OAiPSlw83tqQh5xy+3n+JkcHIj2O71nu6OTwk6Ry0V+6looJ5dh2qgpGE58lTh0UaPdpc1cgs6gWE17cgA1dHiyGhfnjt8dn4d3bJrrkNDUFgx6ELB2vkKjnUyuP4XxyW82A164dq8h1CYLwLO6fnYhxQ4IBAMZOM5749ig67FgulmoMfptZMuB+P+l08oGiOjS0ddi8N4LnN9k9ybOzgmazGU9/fwzX/m+PTIvXR6PC81eOwrY/znVqT2BfUDDoITQZOrArv0ZYK5GOr2w0YGuu2Jx9xbhYl3yiIQjC+WjUKvxj2ThBq+/kuUaZLqnSXD4uBjot/1455U04eW5gmcioIB3GWoLYTjOH7acGPpBCiBhNZmzKEbNjnuw6Ulrfiote3YKvD8j7VyfFh+LIXy/G7dPjnbOxfkDBoIewNbcKRstTeFpMEOLCbE/HP/ndUUlWkOHFJaNtviZBEJ7LiKhAPLUwRVj/Z0u+IHWlNEE6LS4dHSOsv7VhkERaKpY+ABMDZ29hDZosGbLYED+kWYaMPI31J8ox+7WtKGswCH+nUTG8fPUYrLx/ukv1BfYGBYMegtKinmX1bdhxSmz8vSp9sEtMPBEE4dr8bsZwZAzjfX9NZg5/+v6Y3aaLr5NoDq46cg6Gjp7tvHpjjqRUvO1Uld3Fs72B37Lk9yTGmBN3Yx/+8tMJ3PfFQdmkcGKkHrv/PA83TRnqxJ31HwoGPYB2U6fsaXbRaNvT8U9Js4IqhhevpqwgQRB9o1YxvHrNGPio+dvLsbMN+GjXabu819Th4RgSyruINLR1YEPWwAZAxsQGI9Sfb4Gpbm4XfHSJgWE2c11kzjyrRNxqNOHSt7bj871nZH9/+/R4bHpitlsOWVIw6AHszq9Bczufjh8W7o8UG5tUS+tbsVPSf7h0QqzbpLoJgnA+I6IC8egCcbr4H7+dQpEdNPxUKobrMuSDJANBrWKYlSy6NW2jvkGbOHq2HpVNvKRKWIAPJsaHOXlHylFa34rpf9+M7PIm4e+0aob3bs3A81eOcuLObIOCQQ9A6XT8k98eE15rVAwvuPEPOEEQzuHeWQlCn1i7yYw/fX/MLuXXazJicf4jb2d+NUrr2wZ0ndmSYHBrLukN2sJ6SVZwwcgowbLQ3cksqsXc17ehvlWcOI8O9MX2p+a6/bQ0BYNujtnMYUOW+MG10Eah6ZK6VuwpFLOC12UMgY6yggRB9BOtWoXXrh0r8y7++kCx4u8zJNQfM0bwwtEcB/xw8OyAriPNDB4qrieJmQHCcZxcUsZDXEe+yyzBde/uEQY1AWBqQhj2PD0PMSF+TtyZMlAw6OYcLqlHdTOfjg8P8MH4oaE2Xe+Jb48KrzUqhmevSLPpegRBeC+jY4Nx76wEYf3K2hyUNQwsc9cb12aIgyQrD54dUAYyQu8rk5jZJXHOIKynoKpZsPXz91FjhsThxV15a+MpPPndMZmv8M1ThuKbe6dBpfKMMMoz/hVejLRher6N6fgzNS3Yf7pWWN8wKY6yggRB2MSj85OQEMF7mTe3m/DMjycGLBDdE4tGDUKQRe2guLYV+ySfY/2BSsW2Iy0Rz06OhE6rduJubOfZ1Sfx1sY8Yc0APHtFGl66eozzNmUHKBh0czZI+gUvtjEd/8fvxF5BrZrhr1dQryBBELah06rx6rVjhb6+zTmVWH30nOLvcVV6rLBeOcBBkjkp8iESpYNWb+A3hWXOnMlj3xzBp7uLhLVaxfDJnZNw50XDnbcpO0HBoBtTWNWMgio+Ha/TqoS+mYFQVt8me5q+fmKc4CRAEARhC5Piw3Db1GHC+rnVJ4X2FqWQ2tOtPVGGJkP/e/7GDQlBsB8vMVPR2I4cycQo0TdlDW04erYBAN9mNFci5u1u3PHxfvx0RPS89lGr8OPvp2N2ivv+m3qD7vZujLREPDMpEn4+A0/H//F7+QTx8stG2rQ3giAIKU9dkopYS6N9XWsH/vZzlqLXHx0bhNRBvKyWocOMX46V9fsaGrVK1uNGEjP9Q3pPmpYYLgTW7sYN7+2Raff6adVY/9hMjB0S4sRd2RcKBt0Y6S/exTZMEVc3G7AzT2yWvnp8LPUKEgShKHpfDV5eKvZZrT56TlEfYMYYrptou+bgHOobHDDSe5KtyhbO4qb392JvoVglC9RpsPmJ2RgeqXfiruwPBYNuSnVzOw4W1wEAVAyYb0M6/unvj8s8iJ+/inoFCYJQntnJkViSPlhYL//pxIAt5LpjSfpgaNV8c+Lh4nrkV/a/zCsdIsksqhtQudkbaTR0YK9ElmyBGwaDt3ywF7sLxH9DqL8WO/441yOkY/rC64JBxpieMVbCGOMsf+5w9p4GwubsSmHMPWNYKML1vgO6TkOrEZtyxKffxWMHkdsIQRB2Y/nlaUL5sLi2Ff/anNfHGdYTrvfF/FQxCFmZ2X/NwaggnSCWbTJzsuCA6JmtuVXo6ORvSmNigxET7F4B1G0f7ZM5b4X4a7H1yTkI8fdx4q4ch9cFgwBeBDCkz6NcnN8UKhEvX3UC5yW5VAx48SryICYIwn5E6H3x9KWpwvrdbYU4VaHcoMaySeLH+/eHStEhEQm2FulUsbR3jOgZpdqWnMFdnxzA9lNiq1SIHx8IBntJIAh4WTDIGJsA4CEA+5y9F1toM3ZiZ774ATVQSZlmgwlrj4kyAPNTo7zqh58gCOewbGIcJsXzAvkmM4f/++G4YlZ1s5IiERXIV0qqm9sHFMxJS8XbSWKmT4wms6y/csFI9wkGn1p5VFYdC/bTYPMT3pMRPI/XBIOMMRWAdy3LB5y5F1vZkVcFQwf/tDsiSo/hFkHX/vL8zyfRafmQYwBeucazRDQJgnBNVCqGl68eI/T3ZZ6pw4oBDnx0RaNW4RqpI8kArjthWCgCffl2mdL6NuRXNiuyN09l/+laNBlMAIDYED+MjAl08o6s45W12VgpsS8M9NVg0xOzEab3rkAQ8KJgEMDDACYC+C/HcYedvRlbUCIdbzCa8ONhUUNpZnIEIvQ6m/dGEARhDUnRgbhvVqKwfmVtNqqalNEevE4SDG7Oqez3dbVdJGaoVNw7cvODaDA2cCcsR/H+9gK8u71QWOu0Kqx7fKbX3ge9IhhkjMUCeAFABYDlTt6OTXSaOWyWpLQHGgy+tDYHJklZ5tWlY23eG0EQRH94aN4IDAv3BwA0Gkx4cY0y2oMJkXpMHCaWoX+SPPhai7RUTHqDPcNxHDZm235PciTfHyzBS2tzhLVGxbD6oRmIDfF34q6ci1cEgwD+BSAQwJMcxzU4ezO2cKi4DjUtRgBAZKAv0gcggmkymbHigFg6mTI8zCtG5wmCcC10WjVeXCIOra06cg478pQJvJZ10Rzsb9/fbMkQyf7TtWg1mhTZl6eRVdaI0vo2AECQToPJw8OcvKPe2Z1fjSdXiiYLKgZ8dc8UJEe7R2nbXnh8MMgYuwLA1QC2chz3xQDO7/GPM5CWiBeMjIJK1f99vP5bDoySCbvXrqWsIEEQzmFmkn20BxePjYG/xZUpr7JZsEmzlphgP6RYAgRjp1lm10mISO9Jc1OjoFW7blhxuqoZt3+8X9DVZQD+d0sGJg8Pd+a2+oW9YhLX/a4pAGMsAMC/AXQAeNDJ27EZjuNs7hc0m834dM8ZYZ0eF4Jh4QMbQCEIglACqfbgmZpW/Htzvs3X1PtqsHhMjLCWVkOsZaakb3CHRHqEEJEnKFy3RNzQasQV/94laCECwItXj8bCUQNT4/A0PDoYBPA3AEMBvMlx3ICaUTiO6/GPoymoasbp6hYAgL+PGtMTI/o440L+uTlfmEQGKCtIEITz6ao9+N72QhRW2T7BKy0V/3L0HNqM/cs4zpT0DSpVvvYkztW34eS5RgCAVs1k+oyuhMlkxqX/3InmdrHUf9eM4bh5yjAn7mpg2Csm8dhgkDGWDuBRACXgg0K3Ryo0PSspEjqtut/X+HDnaeF1Wkyg1/dJEAThGiybGIcJQ/keaGOnGc+uPmnzDW5SfCjiLQMqTe0mrD9Z3scZcibHh8FHw98m8yqbcc7SG0fwbMwW70lTE8IRqNM6cTc9s+y9PbLv3bzUSPzl8jQn7sj18NhgEMDbANQAngHALDZ0wh/Jcb6Wv3P5MSJbS8Rf7SsWtKAA4KWrSVeQIAjXQKVi+NtVo3G+DXpHXjXWnehf8NYVxhiu6zJI0h/8fNSYHC8OROzMo1KxFHdwHfm/H47hUHG9sE6O1uOD2yY6cUeuiScHg+fzv58BaOrmz3n+Z1kro2lgJyqbDDhSwv9Aq1UM81Kj+n2NNzeeEl4PC/PH+KGhiu2PIAjCVkbHBuPWqWLp7m+/ZKGl3bYp3qUTYoUAc3dBDUpqW/t1vrRvcDuVigUaDR3YWyh6+bpiv+C3mSX4ar/4ABAe4IPVD14ElcqTQ5+BQf9H3IRN2ZU4XzGZOCwUoQH9U0jfmFUuE179y+UjldweQRCEIvxhYQoiLA4QZQ0G/MvGYZKYYD/MTBJ72b6TOE5Yg/TcnfnV6FTINs/d2ZZbJQxjjI4NwmAXkyc7UdqAP38vSsj4qFX45eEZ0PlonLgr18Vjg0GO4+I5jmM9/ZEceqfl7+KdtVdrsDUd/8KabOF1hN4XCwboZ0wQBGFPgv20ePpS8WH1gx2FyK9s6uWMvpEOknx38Gy/fJBTBwUiQs97Hde3duDkObeWqlUMV54ibmg1Ytm7e3D+28wAfHLnJNLT7QWPDQY9iZZ2E3bmi70qC/sZyB0ursOZGrE08oeLkxTbG0EQhNIsnRCLSfGig8hfV9k2TLIgLQoh/vxwQ2l9G/ZIypt9oVIxucQM9Q2io9OMLbmu6TpiNptxxb93olUyOf7nS1MxfUT/1Te8CQoG3YAdeVUwmng5mJToQAwN79+sy/KfTgivA301uMkNx+kJgvAeGOOHSdSWZr/dBTX45VjZgK/nq1FjSXqssF7Zz0ESWd8gWdNh/+laYRgxNsQPaTFBTt6RyL2fH0RxrTg5fMmoaNw3O7GXMwiAgkG34DcbSsRnaloEHSgAuGtmvFLbIgiCsBsjY4Jw+7R4Yf3imiyZTlx/uTZjiPD61xPlaGjrsPrcGZKs0qHiOpv24Ql0bVtyliNXV/6zJV/mkzw8wh/v3DzBiTtyH7w2GJT0D37i7L30hqnTjM05A0/HP/3DceG1j0aFR+ZRiZggCPfgsYuTEBnI9+tVNLbjbYkiQn8ZHRssZLDaTWb8fPSc1edGBemQOojXZO3o5LCvH2VmT6OrE5ar9AsePFOH/7c+V1gH+Kqx+qEZNDlsJfR/ycXJPFOH+lb+CTY6yBdjYoOtPre22Yg9BeKH1rKMIfSLQRCE2xCk0+KZxeIwyUe7ipBbPvBhkusmitnBlf2cKp4lcyPx3r7B7LImlFoEnAN1GkxJCOvjDPvTbDDh1g/3CZ7DKgasvH+ay4pguyIUGbg4XZ/AVCrr0/HLVx0XfjnUjOGZy0hOhiAI9+Kq9MGYMpwPODrNHJ5dfWLAwyRL0mPho+Zve0dL6nGqwvrAkvQGeaT3pLkpUdCqnR9GXPe/3bKBkb9dOQppMdYnTggKBl2arun4/pSIW40mrD8hnrtwVDT8SF+JIAg3gzGGF5aIwyR7C2vx6wCdSUIDfLAgTRTs788gyaT4MPharOkKq1pwtq5/4tWewoZs8f/9AheYIn529UlkS7LFC0ZG4RZJrylhHRQMujCnKppRbFHL1/tqMC0x3OpzX16TjU7L0zMD8NISsp4jCMI9SY4OxG3TRBWEl9Zko02SCeoPUnu6Hw+XoqPTbNV5Oq0ak4d7tzXdufo2nCjlBxK1aoY5KZF9nGFfNmVX4NPdRcI6OtAX792a4bwNuTEUDLowUhPw2cmR8NWorTrPbDbL+mGmJoQjTN8/xxKCIAhX4rEFyQizOC+V1rfh3e0FA7rOrKRIRAfxQynVzUZskQzoWXPuebyxb3CT5J40NSEcQU7syatsNOCBLw4Ja42K4YcHp1Nf/ACh/2sujKxfMM16L+J/bclHu0l82n1lKWUFCYJwb4L9tHhyYYqw/u/WggGVatUqhmsmiIMk32ZaP0gyM1nsG/RGa7rfXGSK2Gw2Y+l/d8Moyeq+dX06YkP6p8FLiFAw6KJUNbXj6Nl6APyH19wU64PBj3aeFl6nxQQhPiJA8f0RBEE4musnxWHUYFEe5pW1OQO6jlRzcEtuJSqbDFadlxIdKEjdNLR14Hip91jTNRo6sFciqePMfsE/fX8cZ+tEYenrJg7B5eMGO20/ngAFgy7KlpxKnB+YyxgWihB/68q8q46UoqFNFER9/spR9tgeQRCEw1GrGJ6TfKatOV4mk8+yloRIvWB312nm8NPhUqvOY4x5rRvJttwqdHTyN6VRg4MQ6ySf3805lbI2qOER/nj92nFO2YsnQcGgi7IhW5qOtz4r+LpEdHNwsA6ThjtfA4ogCEIpJsWH4UpJFuj5n0/CZOUQiJTrMsRBkpWZZ62Wq5H2De7K956+wY3Zzi8RN7Qa8cAXB4W1Vs2w8v5pTtmLp0HBoAti6OiUTapZ+4t3tKROljp/YmGy4nsjCIJwNk8vToWflh+oyylvwtf7i/t9jcVjY+Dvw18jr7IZR0rqrTpv+ghR1eFQcR1ajZ5vTdfRaZYN2vTXCUspbnx/r6wf/s1l6YjQ65yyF0+DgkEXZE9BDdo6eNmEhIgAJETqrTrvr6tOCq8DdRpcI3nyJQiC8BRigv3w4NxEYf2PDadQ32rs1zX0vhpcNiZGWFvrSBIVqENKtGhNt/90bb/e1x05UFSLRgMf9A4O1gl9m47kjd9ykVUm6gkuHjOI+gQVhIJBF0RaIp5vZYm4rL4NR8+Kzcy3k+gmQRAezN0zExAXxvet1bd24I0N/fctlmoO/nzknNXahReNEPsGvaFUvClbzArOHxkNxqx3wlKCE6UN+NfmfGEdoffBv28c79A9eDoUDLoYHMdhs+QXz9oS8V9XnRBea9UMj81PUnxvBEEQroJOq8byy9KE9Rd7zyC7rLFf15gUH4r4cF6OpKndhPUnrXM2mZEklop35vd/gMXd2CwpEc/rRw+7EhhNZtz8wV6Z7/CK+6aRnqDC0P9NF+PkuUaUN/IyB8F+WmQMC+3zHIPRhM054lTb4jEx0GjoW0sQhGezMC0aMyxZOjPHD5P0x7eYMSbLDn5rpT3d5OHh0Fjs8bLLGlHT3N6PXbsXhVXNOF3dAgDw06oxLcF6JywluOezTJlCxtOLU5FoZesUYT0UMbgYUqHpealR0FhhAv7KulyZ9RzJyRAE4Q0wxvDsFWk2+RYvnRALy+nYXVCDktq+haz1vhqkx4UI690DkLdxF6Ql4hlJEdBprXPCUoLVR0uxTSLfMyk+FPfMTOzlDGKgUDDoYmzK6V+/oNlsxrcHxKfZycPDrNYkJAiCcHeSuvEtNnRY71scE+yHmRK5mO+sHCTxlr5B2T0p1XEl4oZWI5789piw9vdR4/PfTXbY+3sbFAy6EOUNBsEEXKNimJXctwn4J7vPCJPHAPDCktF22x9BEIQr0tW3+IMdhf06f5mkVPzdwbMwW2EzJw0Gd3poMNjQ1oEDRXXCep4Dg8GbP9gns5v7780ToPPROOz9vQ0KBl0I6ROYtSbg72wVJ6wSI/VItkgeEARBeAvBflr84WJRV/WdrQVWW8wBvPd7iD//eVta34Y9hX2XfdPjQgSdwrN1bSiu6b9Psquz7VSV4L88dkgwooIco+n3/o4CnDgnDgNdMS4Gs/thyUr0HwoGXYiNWf0rEW/NrUR1s6it9cxlI+2yL4IgCFfnhklxSI7mBwtajZ34x3rrpWZ8NWosSY8V1tYMkvhoVJgicXjyxOzg5mx5D7sjKK1vxd8lntOh/lq8fX26Q97bm6Fg0EVoNZqwS9KEbI2kzAu/ZAuvI/Q+Dk3hEwRBuBIatUomNfPtwRKcPNfQyxlyrs0YIrxed6IcDW0dfZ7jyX2Dpk4ztkqGNxxlQXfje3thsUAGA/Dl3VNIRsYB0P9hF2FnXjWMFpudlOhAxIX593p8fkUTCqqahfXv54yw6/4IgiBcnVnJkZiTwvdacxzw4i/ZVkvNjI4NRloM76zRbjLj56Pn+jxHGgzuLqi2qtfQXThUXI/6Vj4gjg7ydYjryItrslBcK1qq3jVzONIGB9v9fQkKBl2Gjf10HVkuEZn206pxx/RhvRxNEAThHSy/bKQgNbOnsEYm19UX100Us4PW2NOlRAciQs8PrtS1diCrn6LXroy0h31eapTdXUeyyhrw4Y7Twjou1E+W6SXsCwWDLoDZzMlEo+f3kY6vbzViX6Hoh7lsUhyl0QmCIACMiArEzVOGCuuX12YLVZe+WJIeCx+LtuvRknqcqmjq9XiVimF6omeWiqVOWPNT7VsiNpvNuOOjA4LLiJoBX9071a7vScihCMIFOHq2HtUWBfsIvY9MzLQ7nludJfmlYXj6khQ775AgCMJ9eGxBMgJ1vAxJUU0rPttTZNV5oQE+WJAmVmZWWjFIMsMDJWaKa1qRV8m3IflqVLJyuD14/pdsVDaJLi5PLUpFXGjvrVKEslAw6AJIS8RzU6KEEkd3mExmrDku9rLMTY0k7SWCIAgJYQE+eFTiz/7PTXmoazH2coaI1J7ux8Ol6OjsPas4fYRoz3agqBbtJuun5A6qAAAgAElEQVQFr10VaYl4emI4/Hzs5zqSU9aIz3YXCevESD3un0MuI46GgkEXQGr301eJ+J1tBejoFJuUX7iKRKYJgiC6ctu0eMSH89mlRoMJb2/Ks+q8WUmRGGTR06tuNmJLTmWvxw8J9Rfex9BhxqEz9Tbs2jXYnGP9PclW7vh4v6w8/MVd5DLiDCgYdDIlta3IKef7Unw0KsxM6j0d//EuscF2VEwQYkL87Lo/giAId8RHo8LTi0Xt1c/3nkF+ZXMvZ/CoVQxLJ0g1B/seJPEkiZkmQwf2SkS37SlZ9sIvWShvFMvDf1iYQvc0J0HBoJPZlC1Pxwf49lzy3ZRdgbpWUfvqL5fTpBVBEERPLEyLxtQEXhi608zh5bXZfZzBIy0Vb8mt7NPNxJP6BnfmVQvVp7SYIAy2U3CWX9GEj3aKyY2EiAA8OJck0pwFBYNOZlM/0vF//1VUZY8K9MXUxPBejiYIgvBuGGNYflkazquibM6pxI68qt5PAjA8IgCT4kMB8EHkT4dLez1+WmK48B7HztZbJVjtqsjvSfbLCt72kVgeVjHg87upPOxMKBh0Il3T8Qt6+cU7U9MiTHcBwP2zqcGWIAiiL0bHBuPaCaJ+4Iu/ZMPUx1AIIM8Orsw826t4dYi/D0ZbxJHNHLDPCm9jV6TTzMl6JO1VIn5pTRbONYjZ1kcXJCE2hKaHnQkFg05k+ykxHT9qcBBigntOx/911UnhtU6jIpFpgiAIK3lqUQr8LROxuRVNVvUBXjYmRjgnr7IZR0p6HwzxhL7Bo2frUWOZuo7Q+2DckN5lzgZCQVUzPpCUh+PD/fHo/GTF34foHxQMOhFpv2Bvvo+tRpOstHHV+FgSmSYIgrCSqCAdHpBUU97YcAot7aZezwnw1eCyMTHCui9HkoskEjPu2jcoFZqemxIFVS8yZwPltg/343ySlS8PT1H8PYj+QxGFkzB1mrElV/zF6y0YfG1dLs5bXjLG2y0RBEEQ1nP3zASJZEw73tte2Oc50lLxz0fOoc3Ys4bgpPgw+Gj4W2pBVQvKG3ofOnFF+muL2l/e+C0XpfWi9/DD80aQuLSLQMGgkzhUXC9MBkcH+WJ0bM8m4N9KVPCnDA9DoE5r9/0RBEF4En4+avxhoViOfG97ISobew/YJsWHChqCTe0mrD9Z3uOxOq0aE4eFCmt3KxWX1reJMmdqFWYkRSp6/YpGA/6zJV9YDw31w+MXk3uWq0DBoJPYJHsCi+7RBPzr/cVolTyNPnfFKLvvjSAIwhO5ZsIQpA4KBAC0dXTizY2nej2eMSbLDn7bhz2dO/cNSoWmpySEQd+LzNlAuOPj/Tjvl8AAfPw7mh52JSgYdBIbZf2CPafj/7VZfJIaFu6P1JieM4gEQRBEz6hVTCZEveJACU5VNPV6ztIJsTjfOre7oAYlta09HjtdIve1u6Cm1wlkV0OWoFB4ivjbzBJkl4n/n2+aMhSJkXpF34OwDQoGncDp6hYUVLUAAHRaFaYndu86crSkDuck/RVPLKSJK4IgCFuYnRwpOD2ZObl+a3fEBPthpqRk+l0vgyRjYoOFjFp5owFFNT0Hjq5Eq9GE3QWiHI6SFnStRhP+8tMJYR3qr8ULV1GFy9WgYNAJSJ/AZiZFQqft3gT8uZ+zhNdBOg2uHBfb7XEEQRCE9Tx96UiZEPXugt5LusskpeLvDp6F2dx9xk+jVmHK8DBh3dd1XYVd+TUwmnjtxeRoPeLClBvquO/zg2g3ibqO/7slg9QwXBD6jjgBa0rEtc1GHCkWda1umUq6ggRBEEqQNjgIS8eLQtQvr83uMcADgAVpUQjx5wf3SuvbsKcXUelpXUrF7sDmHPGeNC9Vuazgjrwq7MgTA+J5qZGYkkDOWa4IBYMOpqG1AweK6oT13B56M57/5aRg1aNRMTy2gErEBEEQSvHkomT4WqRgTpQ2YvXRcz0e66tRY0m6WJnpbZBEGgzudYO+QbOZw6ZsqcyZMv2CZrMZD311SFjrtCq8c9MERa5NKA8Fgw5m66lKdFqeQNPjQhAVqLvgGJPJjLXHy4T1/NQoQb+KIAiCsJ2YYD/cPXO4sH59fS4MHT3rCF6bIWYS150o79F/eOSgICGLWNNixKmK5m6PcxVOnmtEZVM7ACDEX4vxQ0P7OMM6/vzDcTS0icLef186FjofZSeUCeWgCMPBWPME9r/tBYJNHQA8dyU12xIEQSjN/bMTER7gA4Av/366u6jHY0fHBiPNoubQbjLj5x4yiSoVw7QEaanYtfsGN0lKxHNToqBWwHXkVEUTVkos/8bEBmHJeOp5d2UoGHQgpk4ztkpcR3qa2PpoV5HwOi0mEDEhPXsWEwRBEAMjUKfFowuShPW/t+SjzuLN2x3LJorZwd7s6bpKzLgy0gSFUq4jd358QGhzUqsYPr5zkiLXJewHBYMO5OCZOjQa+LT54GCdIH4qZWtuJWolH0bLL0tz2P4IgiC8jRsnD0VCRAAAoMlgkmm7duWq9Fj4qPnb5tGS+h41CqV9g/sKa4TWIFejotGA46UNAPje9JkKuI68szVfZjn3+IIkROgvbIciXAsKBh2IVOF93siobl1HpJpXEXpfTB/RvQYhQRAEYTtatQp/vCRVWH++twhnalq6PTY0wAcXp4kVnZU9DJIkRuoRGegLAGg0mJB1rlHBHSvHFsk9aVJ8GIL9bLM6bWg14o3fRFeXIaF+eGheUi9nEK6CRweDjLGhjLHHGGM/M8aKGWPtjLEmxthRxtjfGWMxjtzPJskv3vxuxvfLJN6QAHCPpLmZIAiCsA+LRkULvsIdnRxeW5/b47HXSkrFPx4uRUen+YJjGGNdSsWu2Te4UeES8QNfHoJJkgV977YMm69JOAaPDQYZY3EAigC8CeByAHEADAD8AIwF8CcAJxljcx2xnzM1Lciv5KfKdFqVrIxwnud/Pim89lGrKBgkCIJwAIwx/N9lok3dmmNlOFxc1+2xs5IiMSiIL3tWNxtl2TUp0mCwN11CZ2Ho6JT5J9vqOrI7v1rWH3nZmEFIiwm26ZqE4/DYYBDAeVuPNQCuAxDGcVwwAH8AiwGcBhAK4CfG2CB7b0ZaIr4oMeIC1xGTySx7Sls0OppU2gmCIBzEhKGhuGyMWCx6eW12txqBahXD0glSzcHuB0mkNqP7T9d2m0F0JnsKa9BmkdJJiAjAcEvf5EAwm814UKIp6KdV483rx9u8R8JxeHK0UQdgPMdxl3Mc9x3HcXUAwHGckeO4X8EHhAYAQQDus/dmuvYLduWdbQWy9PpfaHCEIAjCofzxkhRo1Xwv94GiOmzIquj2uOsk9nRbcitR2WS44Ji4MH/EWpQgWo2dOHa2/oJjnInUFtXWEvErv+agrlXUXXxxySjSxnUzPPa7xXFcA8dxR3v5eg6AvZalXRsbmttN2FdYK6zndeM68olE32rU4CBEBdH0FUEQhCMZFh4gs/58fX0uTN1k9IZHBGBSPN9j2Gnm8NPh0m6vJ+sbzHedUjHHcdgsqUTZYkFX0WjAhztPC+vkaD2uyYjr5QzCFfHYYNBKzv92qns9ykZ25lXDaPlASYsJQkywXDdwR16VTE7mmcUjQRAEQTieh+aOgN6Xd8rIq2zGD4e6D/Sk2cGVmWe7LSlPH+GafYM55U0418BnMwN1GkyMH7jryN2fHsD5ohZjwId3kKagO+K1wSBjTAPgIsvyhD3fS2oC3l06/pW12cLr8AAfkpMhCIJwEuF6X9w3K0FYv7HhVLc2dZeNiYG/D59HyKtsxpGSC8vA0xLEz/LMM3W92t05EmmJeE5KFLTqgYUCa4+V4XipKJtz0+ShiAv1t3l/hOPx2mAQwIMABgEwA/i0p4MYYz3+sQazmcPmnCph3bVEXNloQFaZKCdz50Xx/fk3EARBEApz18zhiNDzOoHljQZZG895Anw1soGT7hxJBgXrBEFro8mMQz1MKDsauczZwPoFTSYznvpe7MQK0mnwwlVknWpvbI1JesIrg0HG2FgAr1iW/+Y4Lste73W8tAHVzbwJeHiAD8YNCZF9/W+/iG+tVTPcPyvRXlshCIIgrMDfRyOzqXtnSz4aJAMS55GWin8+cg5txgszf1IZsT0uYE1X3dwuZDFVDJidPDDXkad/PI6WdvHf++b16aSA4cZ43XfOIjT9E3i9wYPg9QZ7hOO4Hv9Yg/QJbE5KFFQSE3Cz2Yz1J8uF9YKR0dDQBBZBEITTuWFSHOLD+ZJno8GEd7ZdaFM3KT5UOKap3YR1J8suOEYqMeMKweCWnEqcv31NHBaG0ACffl/jdFUzvpNkQicMDbFZp5CwDltjkp7wqsiDMRYG4DcAwwHkAbiM47gLNQEUpLd+wfe2n0ZHp/gN/OvlJCdDEAThCmjVKjy1SLSp+3hXEc5JPHcBvmTXdZCkK1MTwoTXR0rq0dJussNuracvmTNruOvTTJy/c6kZwwe30dCIu+M1wSBjLBjAegCjARQDWMBxXPciUgpR0WjACUtzLW8CLh8M+UAyjj9yUCBiQuRTxgRBEITzWDxmEMYN4V00jCYz3tp46oJjlk6IxfmCz+6CGpTUtsq+Hq73ReqgQACAyczhQFFt10s4DKPJjO2nxB72BQMIBr/aV4zCatG7+YE5iQjT9z+7SLgWXhEMMsYCAKwFMBFAOfhAsNje7yu1KZqSEIZAnWgCvregRuglBIA/X5oKgiAIwnVgjOFPl4ifzd8dPIu8iibZMTHBfpiZFCk7pivTXMSabt/pGrRY+hqHhvkjMVLfr/ONJjP+9otomxqh98WTi1IU3SPhHDw+GGSM+QH4GcB08LqCCziOy3PEe0v7BbuKer60VhwcCfXXYnaK7SbhBEEQhLJMHxGBWZYhCzMHvLY+94JjlklKxd8dPAuzWd6/5Sp9g5tkQtNR/Z5A/eN3x2DoEEW437mZLOc8BY8OBhljPgB+ADAXQD2AhRzHnez9LGUwdHRiZ57EBFwyvl/dbJBpM902bRgIgiAI1+RPl4jZrw1ZFcjsUupdkBaFEH++8lNa34bdXQK+ycPDhFLyidKGbieT7Q3Hcdgk6WFf0M+Bj9NVzVh1RBTgnpoQhsnDw3s5g3AnPDYYZIypAXwF4BIATQAu5TjuUO9nKcdeqQl4ZADiJSbgL/6SI7zWqBgenpt0wfkEQRCEazBqcDCWpA8W1n//NUc2vemrUWNJeqywXpFZIjs/2E+L0bF876GZ48u1jia/shkltfwATICPGpOHh/Vxhpz7vjgoDo2oGP53i11dXAkH47HBIHh3kWssr7UAfmKMlffw54DSb765B1FPs9mMtSdE+YG5KZEkJ0MQBOHiPLEwBVo1n97LPFOHjZKSKwBcP0ksFa8/WY76VqPs687uG5S2Lc1KjoRPP+47q46U4lRFs7C+d2YCQvxpaMST8OQoRPpv0wGI7uXPwFQ3e4DjOFlvxlxJMPjxriIYTWLPxV+vJMV2giAIVycuzB83TxFbel5bl4NOSW/gyJggjJVMHv90WO5p7Oy+wc2Se1J/NAHNZjOe+fG4sA711+KpRcmK7o1wPh4bDHIct5XjOGbln3gl3/tURTNKLXpUgb4aTIoX0/Hvbi8UXidH68nHkSAIwk14eN4I6H01AHg/4u8PySeHpYMk3xwokZWSJ8WHQmNpHMwpb0KNRE3C3tS1GJF5hu9zZAyYk2J9/uMvq06imZxGPB76jtoBaZPurJRIwQQ8s6gWlU3iB8AfaSSfIAjCbQjX++LeWQnC+s0Np2DoEAOlK9MHQ6flP+9zyptwvLRB+Jq/jwbpcaId6d5Cx+kNbjtVhfNJzPS4EMF3uS/K6tvw9X5RhW1cXDDmkPKFR0LBoB3Y0kO/4ItrsoXXwX5aLEgb5NB9EQRBELZx14zhQjBV1mDAp7uLhK8F6bRYPDpGWK84IB8kmS7pG9xdUA1HIe0X7M8U8X1fHBSCSBUD3rt1otJbI1wECgYVpq7FiINn6gCcT8fzwWB9qxFHLebgAHDzlKFO2R9BEAQxcAJ8NXh0/ghh/Z8t+TKpmGWSQZLVR86hzShmDqdJ+ga7ys/Yi45OM7bmyvUFrWFTdgWOnRUzmzdPGYboIJ3i+yNcAwoGFUaajp8wNBRhFhPwl9Zky8byH1tADbgEQRDuyA2ThyI+nO/3bjSY8M62fOFrU4aHCV9rajfhV4l6xPihIfC1TPGerm5BeYPB7nvNLKpDk4H3Qx4crBOs8XrDbDbjiW+PCutAXw2evzLNbnsknA8Fgwojdx3hn8DMZjN+PnpO+PvZyRH9GusnCIIgXAetWiWzYftkVxHKGvihQcaYLDsoLRXrtGpMGBoqrB2hN7hZ0sM+f2S0Va4jr67LQX2bmO18ZekYGhrxcOi7qyAdnWZsy5WO7/PB4Bf7imGQyMk8ewXJyRAEQbgzi0fHCFIy7SYz3t4oupxeO2EI1JbJ4X2na3G6ukX42tQEid6gA0rFMgu6kX2XiGubjfhgR5GwTh0UiMvHDe75BMIjoGBQQQ6eqUOjJB2fEs2n4/+3tUA4JjEyAMPCA7o9nyAIgnAPVCqGPy5KFdYrD55FYRUvzBwVpMNciXzLtxJHEqn49F47i08XVjWj0BKI+mnVmJbQt33c/V8eRKdFEocBeJ+GRrwCCgYVROo6Mm8kbwKeVdaAc5K+kMepV5AgCMIjmJEUgYtG8AFWp5nDPzacEr4m1Rz8/uBZmDr56tC4uGChb7CoplUoL9sD6T1pRlIEdFp1r8cfPFOH/adFyZsl42MRF05auN4ABYMKsilb0puRyo/vv/iLKCej91VTup0gCMKDeEqSHVxzrAwnLNqCc1OjBAmayqZ2bM2tAsD7GGcMk/QN2lFvUFoinm/FFPEjXx8WXuu0Krx+zVi77ItwPSgYVIgzNS0oqOLT8TqtCtMSw2EwmmRlgKUThjhrewRBEIQdSI8LwaJRonbfa+tzAfBDJtdmiJ/5K6SlYgf0DTa0deBAkRho9iUp8+W+M4JzFgA8felIaGjQ0Wug77RCyNLxI/h0/D82nBJkZhiDrL+EIAiC8AyeXJgCy7wItp+qEpIAyyaKweDmnEpUNvEtQ1OlfYN2mijekVcFk+UGNHZIMKJ60Qg0mcyyKlZMsA63T4+3y74I14SCQYWQ9QtaSsRSSYGJw0Kh12kcvi+CIAjCviRFB+Lq8WLg99q6HHAch4RIPSZbvOk7zRy+P1gKgA/OztvWnalpxbl65fsGZVPEfWQF/7r6BNoktnpvXZ+u+H4I14aCQQVobpeXg+elRmFTdoUwWQwA/7d4pDO2RhAEQTiAxxYkwcfiQ3+ouF4IxqSagyszS8BxHHw1akwcFib8vdJTxZ1mDlukMmepPVvQVTYa8I0kcTFhaAimWDF1THgWFAwqwM68KnR08un4tJggDArW4XVL3wgARAX6YrxEaJQgCILwLOLC/HGTxGb09fW56DRzWDxmEPS+fFWosLoFB4p4u9KpCfYLBg8X16HeYpEXHeSL0bFBPR770FeHZP7D/70lQ9G9EO4BBYMKIJvYGhmFikYDcsqbhL+7a8ZwZ2yLIAiCcCAPzRsBfx9eviW3ogmrj5bC30eDKyQqEufbh6Ti03sVnije2KVE3JPryMEzddhvCU4B4JqMIeQ/7KVQMGgj5i7p+HmpUXhxTZaw1qoZ7plJwSBBEISnE6H3lT38v7HhFIwmM26QlIrXHi9Dk6EDY4eEwM+i+1dc2yqb5LUVqQXdvF5KxA93kZJ5ZckYxfZAuBcUDNrIsdIGVDcbAQDhAT4YMzgI60+Iv4gLRkaTpyNBEISXcM+sBIT4awEAJbVtWHGgGGOHBCN1EO9I1dbRiZ+PlsFHo8LEeLF9aK9CEjMlta04VcE7ofhqVJgxIqLb477cd0Y2uEJSMt4NfedtZLNEaHpuahQ+31cMY6foQ7z8chocIQiC8BaCdFr8fk6isP7n5ny0dXTKHEnOaw7KS8XKBINS84PpieHw87nQdYSkZIiuUDBoI5tz5Qrv720vFNYjovSIDSErH4IgCG/itmnxGGTpvatqasfHu4pw9fhYYdr4aEk9csobZcHgHqWCQZktavcl4r+skkvJvH0DScl4OxQM2kB5gwEnShsB8L2BkXoflEl8iJ9YSD7EBEEQ3oZOq8Yj85OE9bvbCqBiDBdLnEpWHCjB2CHBQt/g2bo2lNS22vS+ze0mmb1ddxZ0lY0GmRtKxrAQTB5OUjLeDgWDNiAdHJkyPBz/T2JSHuirwaWjY5yxLYIgCMLJXDdxCOLD+cpQo8GE/20vwPWSUvGPh0th5jhZ3+C+07ZNFe/MqxLalEbGBGFwiN8FxzzYRUrmnZtJSoagYNAmpJIyM0aEY7/kF/maDPIhJgiC8Fa0ahX+sDBFWH+86zSSo/WItQRo9a0d2JBVgWmJyvkUy2TOuskKZhbVCjqHAHAtSckQFigYHCCGjk7syq8W1gXVLTIf4iclHwIEQRCE93H5mBikxfCCz4YOM/6zpQDXSfyKVxwoUWyIhIO8WjV/5IXB4CPfHBFe67QqvExSMoQFCgYHyJ7CGqEBNyEyAOtPlAtfm0Q+xARBEF6PSsXw1CViYuDr/cWYnhiO8xrQO/OrEeavFYSqS+sH3jdo1McIMmcReh+MGxIi+/oXe4pISoboEfpJGCCbJen4xEi9zId4+eVpztgSQRAE4WLMSY7E5Hjees5k5vDVvmJB+4/jgB8Pn8OkeNGabqBTxa0hCcLruSlRUKlE1xGTyYyX1uYIa5KSIf5/e/caHFddxnH8+yRpmt7SSem9tA2lFGhruUiLUpQCtqjcZRAUUZkiKioqFwfE6gvGC4ijgo46jpcXXl6ggiiOWlMaBqFoEWEsUAqhhV5I03sbmqQkjy/Oye7JZjfZTXf3tGd/n5nMnst/M0/bTc+Tc/7/58mkZHAI3J1VkeX7L2zbm9qeVD+cBRm/kYmISGUyM74UuTv4x2e38q4T0oWgf/f0Zs487vD7FB9sSNc2zHxErFIyMhglg0PwUuuBVOug0cNr2Lw7XU7merWeExGRiDMax6UWdLgH3UYawi4lW/YcTJWXAXiqZRfuXtD3f6t2DF2jgrI1tdVVnH3ChNS5zFIyZ8xsUCkZ6UfJ4BA0Rfo+1kfmBg6rNpYvVjIoIiJ93XrBiam5gqvWt3FWZBXxvzbtYlSfeYOF9SmOPiI+c9Y4Rg9PX5f6l5I5fYh/AkkyJYNDEJ0v2LqvM7W9dK76EIuISH8nT6nnklOmpvY3RRaKND2/nVNnpKcXFfqouM8j4khJmWylZCaqlIxkocylQLvau/jPa+kfru7I7fwVF2rhiIiIZHfz0jnUhAs7/rdlH7PGjwKgq7uHupr0o+JCksGDXd10jJ2R2j8/0oLupt8+k9oeMaxapWQkJyWDBWp+aXvqlvuw6vRqrTmTRjMlS7V3ERERgJnHjOKqhekuJNFFHRu2H0htP9myM+95g/98eQdeFcw/nDNpNNPHBV1PfvXkRrZG2qPe8f6TVEpGctIno0DRCu+HutM/rCoyLSIig7np/BOoGxZcerft7aC2Oth+bdebqYUk2/Z28Fqe9QabIpUtzjspuCuYrZTMR9/ZWIzwJaGUDBbgUHcPzS+19Ts+pq6GZfMmxxCRiIgcTSbV1/Hxs9ILDWsjd+t6VxhDfo+KgzJn6QWNvSVl7swoJXP/h047rJgl+ZQMFmDtxt3sjxSX7nWl+hCLiEiePn3O8YwJK1Ec6ExfU3o7iEB+fYrXbd2XWsRYdeggp89oYPu+Dh7IKCVzRqSotUg2SgYLEO372KvKgpIBIiIi+Rg7chifOie9Arg6XFTS1d2TOrYmj3qD0WlLI/a0UF1l3PhrlZKRwikZLEDTC639ji1sHMfIWvUhFhGR/F23uJHxo4cD0N2TTvqqw2KEb+zrYNPOgecNRmvejtzzCms37mLtpnS1iyvPmK5SMpIXJYN52rijnVfa2vsdX6E+xCIiUqCRtTV87rzZ/Y5Hy5UN1Kd4+74OntsctkLt6aZuz0Y+l1FK5uuXzi9ewJJoSgbzFO1F3GvK2DrmTxsbQzQiInK0u3rRdKYNUJJsoEUk0WlLdfs30z7+ZLZFSsncqVIyUgB9UvKULRm8/my1nhMRkaEZXlPNF95zQs7zawaoN/iPyHzBut0t7J6xJLU/dWwdH1EpGSmAksE87O84xFOv9v0Nrba6iusWN8YTkIiIJMLlp03j+Amjsp5r3dfJqzv6T0/qONTN4xt2pPY7x0zGq9Nlae5TKRkpkJLBPDy+YUefAtMAy+apD7GIiByemuoqbhmgacGall1Zju1M1RGcMW4kB8el37+wUaVkpHDKZvLQlOURsRaOiIhIMbx33mTmT6vPei7bvMFoSZmOQ91gwaW8yuCHH1YpGSmcksFB9PQ4qzPqC544aQyTtFxfRESKoKrKcrY0zexTHHQdSV+Ttu/vTG2rlIwMlZLBQTy3ZW+fqvAAtyybE1M0IiKSROfMmcCiLI932/Z30hKZN7i+dT9b9hwEwCLjrLtLpWRkyJQMDmJVRqHpevUhFhGRIjOznN2soo+Ko4+IozPZGzatVikZGbKK+OSY2WQz+76ZvWJmHWbWamZ/MrPzB3vvyuf7JoMfXDi9ZHFKZWhubqa5uTnuMKTI9O8qh2vRceNYcuKEfseji0iydcJ6a+92nnngvpLGJsmW+GTQzBYA/wNuAmYBncB44CJgpZndPtD7X3hjf/p7ATcv1SNiEREpjWxzB594eQfuzs4DnTzz+p5+59sevqccoUmCJToZNLMRwMPAMcAzwHx3Hws0AN8hyO++YWbL8vl+75h1jPoQi4hIycyfNpYL3zalz7Gd7V280tbOo+vbyKxBvbCxga6tL5YxQkmiRCeDwCeBmcAB4GJ3Xwfg7vvc/VbgIYKE8Jv5fLMVF51cqjhFREQA+OLSOX0WhwCsaWAMteEAAAWSSURBVNnBqhf7PiJWKRkplqQng9eEr79x9y1Zzn87fD3dzHJX/QQmjK5l7lT1IRYRkdKaPXE0l5wytc+xPz+3jUdfbOtz7CqVkpEiSWwyaGZjgLeHu3/LMWwNsDfcHnAxyY3nzi5SZCIiIgO77b1970+sadmV6joCUFdTxV0qJSNFkthkEDiZdBmmddkGuHsPsD7czdlSpMrgY++cWdzoREREcji2YSQXzJ2U8/yKi+eqlIwUjXnmbNSEMLNLCeYEAtS7+/4c4x4ELgP+4O5XRI4n8y9GREREEsvdM6ecDirJv1aMimwfHGDcm+Hr6BLGIiIiInJEUp2UHIaSWYuIiIgcbZJ8Z7A9sj1igHEjw9cDJYxFRERE5IiU5GRwa2R7as5R6XPbShiLiIiIyBEpycngi6T7eM/LNsDMqoDe9fvPlyMoERERkSNJYpPBcPXw2nB3aY5hZwK9laSbzOyXZuZ5fi0p9Z9Bjm5mtsDMusLPy/JBxt4bjms1s2PKFaMMXeT/i9UDjJloZs+G43aZ2cIyhigJpuuVQH7/D+Ujsclg6Dfh6zVmNiXL+VvD16fdfX3k+CGgdZCvrpJELInh7s8Bd4e79+b4DGJmi4AvhLufdfed5YhPSiv8914NLADagHPd/d+xBiVJpOuVHLakryb+CcFFdibwZzO71t2fD7uTrAA+EI77csb7nnD3JeULUxLsLoLP2Vzgh6Q/cwCYWS3wc6AaeMjdHyh7hFJ0ZjYdWAXMJpiPfL67vxBvVJJQul7JYUv0nUF3PwhcCuwETgfWmdleYA9wG8Gcwjvc/e/xRSlJ5u5dwHKgB7jczK7MGPIVgjmte4AbyxyelICZHQc8RpAIvg6co0RQRI5kiU4GAdz9WWA+cB/QAgwnSA4fAZa6+7diDE8qgLuvAe4Pd39gZuMAzOwU4Pbw+M3urhXtRzkzm0OQCDYCrwLvdvcNsQYlIjKIxCeDAO7+hrt/3t2Pd/c6d5/o7he5e1PcsUnFuJMgOZgIfM/MaggeDw8DVrr7L+IMTg6fmc0DmoFjgZcIEsGNsQYlIpKHikgGReLm7u3AJ8Lda4HfE0xdaAduiCsuKY7wLu9qYDKwjuDR8OZYgxIRyZOSQZEyCe9E/yzcvSR8vUN3j456jcCjwHjgv8ASd38j1ohERAqgZFCkvH4U2d5AsMJYjm4zgYZw+zPuviPOYERECqVkUKRMzMyAeyOHZgOLYwpHiqcF6E0AHzCz2XEGIyJSKCWDIuVzA7AE6CSYX2bAT81seIwxyeF7HXgPsJug13mTmc2MNyQRkfwpGRQpAzM7Frgn3L0LuJogeTgR+GpccUlxhCWslgF7gRnAKjObFm9UIiL5UTIoUh4/BuqBZ4G73b0VuCU8d5uZLYgtMikKd18LvA84AMwiuEM4Kd6oREQGp2RQpMTM7BrgQqAbWO7ubwGEtQWbCGoN/tTM9PN4lHP3J4GLgDcJ7vo2mdn4eKMSERmYLj4iJWRmE4DvhbvfcfenM4bcQJA4LAI+X87YpDTcvRm4jGBu6DxgpZk1DPwuEZH4KBkUKa37CerPbQC+lnnS3Vsix+8ys8ayRSYl4+4rgSuALuBU4K9mVh9vVCKSYMPMbPwgX8NyvVnJoEiJmNklwFWAA9e7e0eOod8FngZGEcwtlARw90cIFgq9RXDn9y9mNireqEQkoc4C2gb5ylnKTMmgSAmY2VjSBaZ/4u6P5Rrr7t3AcoKk4QIzu7YMIUoZuPuDwEcI5osuBv5kZiPijUpEpC9z97hjEBEREZGY6M6giIiISAVTMigiIiJSwZQMioiIiFQwJYMiIiIiFUzJoIiIiEgFUzIoIiIiUsGUDIqIiIhUMCWDIiIiIhVMyaCIiIhIBVMyKCIiIlLBlAyKiIiIVLD/A1yhXdtx+vACAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 720x720 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "figure(figsize=(10,10))\n",
    "set_fig_properties([gca()])\n",
    "vlines(sym_points, ymin=0, ymax=17)\n",
    "plot(linear_path, nu, color='C0',lw=3)\n",
    "xlim([0, max(linear_path)])\n",
    "gca().set_xticks(sym_points)\n",
    "gca().set_xticklabels([r'$\\Gamma$','X', 'K', r'$\\Gamma$', 'L'])\n",
    "ylim([0, 17])\n",
    "ylabel(r'$\\nu$ (THz)')\n",
    "show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Phonon dispersion of silicon crystal described by the mini-Tersoff potential."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "- The above figure shows the phonon dispersion of silicon crystal described by the mini-Tersoff potential [[Fan 2020]](https://doi.org/10.1088/1361-648X/ab5c5f)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 4. References\n",
    "[Fan 2020] Zheyong Fan, Yanzhou Wang, Xiaokun Gu, Ping Qian, Yanjing Su, and Tapio Ala-Nissila, [A minimal Tersoff potential for diamond silicon with improved descriptions of elastic and phonon transport properties](https://doi.org/10.1088/1361-648X/ab5c5f), J. Phys.: Condens. Matter **32** 135901 (2020)."
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}
